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ABSTRACT 


The  servo  actuating  and  corresponding  energy  distribution  subsystem  com¬ 
prise  approximately  80%of  an  aircraft's  flight  control  system  weight.  Con¬ 
sequently,  whenever  subsystem  tradeoff  studies  are  conducted  it  is  desirable 
to  select  the  optimum  design  with  respect  to  weight  and  other  similar  param¬ 
eters.  This  study  investigates  and  develops  such  an  optimal  design  process. 
A  sample  problem  was  selected  an  an  optima1  technique  iormulated  and 
applied  to  the  problem.  The  sample  problem  was  a  fixed-configuration 
hydraulic  actuation  and  power  system.  The  study  objectives  were  to  optimize 
weight,  dollar  cost,  size,  dynamic  performance,  and  reliability  as  a  function 
of  the  system's  independent  design  parameters.  The  parameters  included 
pressure,  actuator  area,  actuator  torque  arm,  and  plumbing  tube  sizes. 
Parameter  optimization  was  accomplished  by  fixed  grid  and  random  search¬ 
ing  techniques.  Within  the  framework  of  parameter  optimization,  a  design 
philosophy  was  formulated  which  allowed  dissimilar  terms  (e.  g.  ,  weight  in 
pounds  and  dollar  cost  in  dollars)  to  be  combined  to  form  a  total  performance 
criterion  for  the  system.  When  the  optimization  technique  was  applied  to  the 
sample  problem,  the  performance  criterion  showed  little  variation  as  a  func¬ 
tion  of  the  parameters  being  optimized.  All  of  the  cost  functions  had  large 
nominal  values  and  only  slight  variationr  about  this  nominal.  To  realize  the 
potential  of  the  design  technique  developed  in  this  study,  different  design 
concepts  and  possible  variations  of  each  one  should  be  considered,  to  reach  a 
more  meaningful  optimum  design. 
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SECTION  I 


INTRODUCTION 


1.  BACKGROUND 

For  each  advanced  aerospace  vehicle,  tradeoff  studies  are  performed  to 
develop  the  best  design  for  the  established  requirements.  So  that  the  trade¬ 
off  studies  are  meaningful,  the  subsystem  and  system  designer  must  have 
realistic  information,  particularly  on  weight  and  performance.  In  the  area 
of  sei  vo-actuating  subsystem  design,  this  information  often  is  based  solely 
on  an  individual's  educated  guess  rather  than  on  an  established  technical 
method.  This  policy  allows  the  bias  from  the  designer's  background  and 
experience  to  influence  the  system  design;  also,  often  it  is  difficult  to  deter¬ 
mine  whether  the  selected  design  is  optimum  with  respect  to  those  factors 
that  are  used  to  judge  a  particular  design,  such  as  weight,  performanr  e, 
size,  cost,  and  reliability. 

The  purpose  of  this  study  was  to  investigate  the  selection  of  parameters  for 
an  optimum,  actuation- system  design  for  a  fixed -configuration  system.  This 
optimization  procedure  is  commonly  called  "parameter  optimization.  "  This 
is  only  one  small  portion  of  the  optimal  design  problem.  For  a  truly  optimal 
design,  it  would  be  necessary  to  consider  different  design  concepts.  Within 
the  framework  cf  parameter  optimization,  a  design  philosophy  was  formu  ¬ 
lated  that  allowed  dissimilar  terms  (e.  g.  ,  weight  in  pounds  and  dollar  cost 
in  dollars)  to  be  combined  to  form  a  total  performance  criterion  for  the 
system.  The  results  of  this  work  were  applied  to  a  sample  problem. 

This  study  concentrates  on  the  muscle  portion  of  the  servo -actuating  subsys¬ 
tem;  that  is,  the  energy  source  and  the  valve  actuator.  The  reason  for 
deciding  on  this  subsystem  was  that  the  majority  of  the  flight-control  system 
weight  is  associated  with  the  actuation  and  power  supply  system.  For 
example,  the  estimated  weight  breakdown  for  Douglas's  proposed  C-5A 
flight  control  system  is  as  follows; 


Pilot 

Con¬ 

trols 

Auto¬ 

matic 

Flight 

Control 

System 

Actuation 

Systems 

Flight  Control 
Hydraulic  System 

3% 

16% 

40% 

41% 

With  regard  to  the  titles  the  above  chart,  PILOT  CONTROLS  and  AUTO¬ 
MATIC  FLIGHT  CONTROl.  J'VSTEM  are  self-explanatory;  however, 
ACTUATION  SYSTEMS  and  FLIGHT  CONTROL  HYDRAULIC  SYSTEM  need 
additional  explanation  Actuation  Systems  rcfci  iu  all  components  that  per¬ 
tain  directly  to  the  aileron,  elevator,  and  rudder  actuation  systems.  Flight 
Control  Hydraulic  System  refers  to  all  components  primarily  in  the  hydraulic 
power  supply  system  that  do  not  pertain  directly  to  any  individual  actuation 
system. 
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2.  APPROACH 


The  following  paragraphs  describe  the  organization  of  the  plan  used  to 
develop  the  optimal  design  process. 

a.  Determination  of  the  Present  State-of-the-Art 

A  survey  was  conducted  to  determine  the  status  of  actuation  system  design  in 
the  aerospace  industry.  The  actuation  systems  were  restricted  to  those 
having  torque  outputs  and  included  hydraulic,,  pneumatic,  and  mechanical 
systems  used  in  the  industry. 

b.  Selection  of  a  Sample  Design  Problem 

Based  on  the  data  obtained  from  the  survey,  $n  actuation  design  problem  was 
chosen  for  use  in  the  development  and  aoplication  of  the  optimization 
technique. 

c.  Generation  of  a  Mathematical  Model 

Because  the  optimal  technique  is  an  analytical  process,  it  was  necessary  to 
generate  mathematical  relationships  of  the  cample  problem  correlating  its 
design  parameters  with  the  system  outputs  being  optimized. 

d.  Formulation  of  a  Specific  Optimization  Technique 

The  parameter  optimization  process  best  suited  to  this  problem  had  to  be 
selected  and  a  method  formulated  for  combining  the  various  system  outputs 
so  that  meaningful  results  could  be  obtained. 

e.  Integration  of  the  Optimization  Technique  With  the  Design  Problem 

A  theoretical  system  design  utilizing  the  sample  problem  and  the  desired 
optimization  technique  was  demonstrated  with  the  aid  of  a  computer.  In  this 
way,  general  design  guidelines  were  established  and  solutions  to  major 
problem  areas  could  be  emphasized. 

f .  Evaluation  of  the  Derived  Design  Technique 

The  usefulness  of  the  new  design  technique  was  evaluated  by  comparison  with 
current  design  methods  and  with  any  related  efforts  outside  of  Douglas.  A 
review  council,  comprised  of  experienced  Douglas  design  engineers,  was 
established  to  provide  counseling  during  the  study  and  evaluation  of  the 
results  upon  completion  of  the  study. 
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SECTION  II 


DESIGN  PROBLEM 


1.  STATE-OF-ART  SURVEY 

The  actuation  system  statg  of  the  art  was  established  by  a  literature  search. 
Its  purpose  was  to  determine  the  state  of  the  art  of  servo  actuation  techniques 
in  contemporary  aerospace  flight -control  systems.  The  ground  rules,  made 
to  regulate  the  extent  to  which  the  various  aircraft,  missile,  and  booster 
actuation  systems  were  investigated,  are  as  follows: 

a.  Only  those  systems  that  had  torque  outputs  were  considered. 

b.  The  choice  of  an  actuating  system  often  depends  on  its  required 
power  supply;  therefore,  the  po  .er  supplies  also  were  included. 

c.  Only  the  valve  actuator  and  load  were  considered  in  the  servo 
actuation  portion  of  the  system.  Amplifiers,  feedback  tranducers, 
and  otlier  electrical  components  fall  logically  into  another  category 
and  were  not  included.  The  power  distribution  and  usage  portions 
were  the  elements  included  in  the  survey  (Figure  1). 


Figure  1.  General  Flight  Control  Actuation  System 
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The  type  of  systems  investigated  were  divided  into  three  basic  design  concepts: 
(1)  electrohydraulic,  (2)  electropneumatic,  and  (1)  electromechanical. 

Because  of  the  amount  of  information  collected,  it  is  impossible  to  include 
all  of  it  here.  Essentially,  the  information  consisted  of  a  general  comparison^ 
a  review  of  current  and  advanced  technology,  and  a  discussion  of  problem 
areas  for  each  concept.  With  regard  to  hydraulic  systems,  advanced  tech¬ 
nology  is  primarily  concerned  with  reliability  improvement,  efficiency 
improvement,  and  compatibility  with  high-temperature  environments.  Its 
principal  problem  areas  are  associated  with  temperature  sensitivity,  contam¬ 
ination  sensitivity,  storage,  serviceability,  and  multiple-energy  conversions. 
With  regard  to  pneumatic  systems,  compressibility  is  the  major  problem 
area.  The  advanced  technology  is  primarily  concerned  with  the  reduction  of 
actuator  volumes  and  the  use  of  various  types  of  feedback  compensation.  With 
mechanical  systems,  clutching  is  the  major  problem  area,  and  as  a  result, 
these  systems  are  limited  to  those  having  relatively  small  loads.  Therefore, 
advanced  technology  is  mainly  concerned  with  the  development  of  clutching 
arrangements  capable  of  handling  larger  loads. 

The  results  indicated  that  hydraulics  is  the  most  universal  design  concept 
and  the  most  widely  used  type  of  actuation  system  in  the  aerospace  industry. 
Hydraulics  has  the  unique  capability  of  combining  the  flexibility  of  pneumatics 
with  the  stiffness  of  a  mechanical  device,  thus  making  it  adaptable  to  most 
applications.  The  extensive  use  of  hydraulic  actuation  systems  and  the  large 
quantity  of  design  information  available  on  hydraulics  led  to  the  recommen¬ 
dation  that  a  hydraulic  concept  be  used  for  the  optimization  study  design  prob¬ 
lem.  Because  of  the  greater  familiarity  with  hydraulics  in  the  field  of 
actuation,  fewer  problems  were  anticipated  in  the  application  of  new  design 
techniques.  Furthermore,  a  design  study  based  upon  the  hydraulic  concept 
will  yield  a  more  meaningful  evaluation. 

2.  SAMPLE  PROBLEM 

The  flight  control  systems  of  piloted  aircraft  are  undergoing  an  evolution. 
Historically,  the  reliable  primary  flight  control  system  consisted  of  a 
mechanical  link  between  the  pilot  and  the  control  surface.  With  the  advent 
of  supersonic  and  hypersonic  aircraft,  however,  it  is  becoming  more  difficult 
to  cope  with  the  problems  associated  with  a  mechanical  linkage  system,  such 
as  backlash,  friction,  inertia,  elasticity,  and  other  nonlinearities.  As  a 
result,  considerable  effort  is  being  expended  in  the  development  of  alternate 
control  schemes  with  adequate  reliability,  that  is,  fly-by -wire  systems.  In 
this  type  of  control  scheme,  redundancy  is  used  to  obtain  the  reliability  pre¬ 
viously  associated  with  a  mechanical  linkage  system.  An  example,  that 
emphasizes  the  problems  associated  with  a  mechanical  linkage  system  is  the 
flight -control  system  proposed  by  Douglas  for  the  C-5A.  This  system 
includes  a  fly-by-wire  mode  of  operation  so  that  it  can  meet  the  teiiain- 
to.Uowing  requirements.  These  requirements  are  so  severe  chat  it  is  not 
possible  for  the  mechanical  linkage  system,  by  itself,  to  provide  the 
necessary  performance. 
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The  newly  designed  subsystems  and  components  for  a  fly-by -wire  system  will 
have  many  stringent  requirements,  such  as  high  reliability,  low  cost,  mini¬ 
mum  weight,  and  dimensional  limitations,  and  must  be  able  to  operate  nor¬ 
mally  over  a  wide  environmental  range.  These  factors,  plus  the  fact  that  the 
actuation  and  power  supply  systems  dominate  the  flight -control  system  weight, 
resulted  in  the  decision  to  use  an  aircraft  system  for  the  design  problem. 

Also,  because  Douglas  had  recently  proposed  a  flight-control  system  for  the 
C-5A,  it  was  decided  to  use  the  C  -5A  hydraulic-system  design  requirements 
as  a  guide  in  establishing  the  problem.  Generation  of  an  optimal  design  tech¬ 
nique  was  considered  more  important  than  optimizing  a  particular  system,  so 
it  was  further  decided  to  select  only  a  portion  of  the  flight-control  actuation 
system  and  to  simplify  it  to  only  the  fundamental  characteristics.  After  an 
optimal  technique  and  corresponding  procedural  guidelines  have  been  estab¬ 
lished,  recommended  future  work  will  consider  optimization  of  a  a  actual  sys¬ 
tem.  In  line  with  these  decisions,  a  simplified  version  of  the  C-5A  elevator 
actuation  system  was  chosen  for  the  design  problem.  A  schematic  of  this 
simplified  system  is  shown  in  Figure  2.  Its  essential  features  include  the 
following: 

a.  A  tandem  valve  actuator  combination. 

b.  Dual  control  surfaces,  thus  requiring  four  separate  valve  actuators. 

c.  One  hydraulic  power  supply  consisting  of  a  pump,  reservoir, 
accumulator,  and  a  considerable  length  of  plumbing. 

d.  Mechanical  linkage  between  the  pilot  and  control  surface  as  a  backup 
for  electrical  input  signals. 

The  vehicle  performance  requirements  that  apply  to  the  elevator  system  and 
used  in  this  study  are  the  following: 

a.  A +25°  to  -15°  maximum  control-surface  deflection. 

b.  A  maximum  control  surface  rate  of  40°  per  sec. 

c.  A  maximum  hinge -moment  capability  of  48,  000  in.  lb  per  actuator 
per  piston. 

d.  Design  considerations  at  normal  operating  fluid  temperatures 
include  a  maximum  rate  at  zero  hinge  moment  and  oue-half 
maximum  rate  at  three-quarters  maximum  hinge  moment. 

O 

e.  A  design  consideration  at  -40  F  fluid  temperature  of  one-quarter 
maximum  rate  at  zero  hinge  moment. 

f.  A  bandwidth  r>f  IQ  rad/ sec  ar.d  an  over  damped  step  ic^uuse  wilh  a 
0.  1  sec  rise  time. 

In  this  design  problem,  the  requirements  will  remain  fixed  and  will  constitute 
the  minimum  performance  that  the  actuation  system  must  meet,  regardless 
of  its  design. 


The  actuation  system  dependent  variables  that  will  be  optimized  and  that  are 
called  vehicle  constraints  are  the  following: 

a.  Weight. 

b.  Dollar  cost. 

c.  Size. 

d.  Reliability. 

e.  Dynamic  performance. 

These  constraints  were  selected  arbitrarily.  Other  factors  could  be  added  to 
this  list;  however,  these  five  seemed  to  be  the  most  pertinent  at  the  time  the 
study  originated. 

The  following  actuation  system  design  parameters  that  will  be  considered  as 
independent  variables  in  the  optimization  process  are: 

a.  System  pressure  having  a  range  of  values  from  2,  000  psi  to 
6,  000  psi. 

b.  Actuator  piston  area  having  a  range  of  values  from  nominal  (that 
area  required  to  meet  the  performance  requirements)  to  a  +2  sq  inch. 

c.  The  length  of  the  torque  arm  between  the  actuator  and  the  load 
having  a  range  of  values  from  4  in.  to  8  in. 

O  o 

d.  Fluid  temperature  having  a  range  of  values  from  -40  F  to  +200  F. 

e.  The  plumbing  tube  sizes  have  a  range  from  0.  5  to  1.0  in.  OD„ 

3.  MATHEMATICAL  MODEL 

For  application  of  an  analytical  technique  to  the  design  of  hardware,  mathe¬ 
matical  relationships  must  be  formulated  that  relate  the  hardware  design 
parameters  to  the  end  product.  In  this  case  the  enc  product  is  a  theoretical 
design  having  an  optimum  combination  of  weight,  cost,  size,  reliability,  and 
dynamic  performance.  For  example,  the  weight  of  each  component  being 
considered  will  be  related  to  the  applicable  system  design  parameters.  A 
combination  of  all  the  necessary  relationships  produces  the  mathematical 
model.  When  combined,  the  relationships  which  are  ^dependent  of  time,  are 
termed  the  steady-state  model.  Those  which  depend  on  time  are  called  the 
dynamic  model.  The  boundarv  mnriitinna  «r  assumptions  used  to  guide  the 
formulation  of  these  mathematical  models  include: 

a.  The  hydraulic  system  state  of  the  art  is  based  upon  a  nominal  supply 
pressure  of  3,  000  psi.  For  pressure  greater  than  3,  000  psi,  a 
penalty  factor  is  added  so  that  approximately  the  same  factor  of 
safety  is  retained. 
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b.  The  same  component  configuration  will  be  used,  regardless  of 
variation  in  the  independent  design  parameters. 

c  The  effect  of  the  pump  power  requirement  upon  the  vehicle  will  not 
be  considered. 

d.  The  valve  actuator  supporting  structure  will  be  considered  a  rigid 
body. 

e.  The  torque  arm  between  the  actuator  and  the  load  will  be  considered 
a  rigid  link  with  its  length  being  the  only  item  of  interest. 

f.  The  effect  of  development  time,  experience,  quantity,  expected  life, 
and  detail  design  will  not  be  considered  on  any  of  the  vehicle  constraints. 

g.  A  size  factor  will  be  used  instead  of  volume  because  the  volume  of 
some  components  is  not  too  meaningful.  This  characteristic  has  to 
be  judged  on  the  basis  of  the  type  of  vehicle  being  used;  for  example, 
transport  versus  fighter  vehicle. 

h.  The  electrical  portion  of  the  flight -control  system  is  not  being 
considered. 

i.  The  dynamic  response  of  a  hydraulic  actuation  system  is  not  the 
same  for  different  input  amplitudes.  As  a  result,  there  could  possi¬ 
bly  be  an  optimum  dynamic  performance  associated  with  each  input 
amplitude  and  so  the  decision  was  made  at  the  beginning  of 

the  study  to  simplify  this  requirement  and  limit  investigations  to 
10%  inputs. 

a.  Steady-State  Model 

The  steady-state  model  is  separated  into  sections,  each  of  which  represents 
an  individual  component  or  subassembly.  The  required  relationships  are 
developed  and  the  corresponding  output  equations  for  each  section  are  sum¬ 
marized  in  Appendix  I.  Each  of  the  various  component  and  subassembly  con¬ 
figurations  will  be  described  briefly.  Obviously,  more  or  less  sophisticated 
relationships  could  have  been  used  and  different  assumptions  made.  However, 
the  fact  that  they  have  been  considered  and  used  is  of  prime  importance 
because  they  provide  the  ground  work  upon  which  the  optimization  structure 
must  be  built. 

( 1 )  Pump 

The  pump  is  a  variable-delivery,  variable-angle  unit  with  a  rated  speed  of 
3,  6C0  rpm.  A  photograph  of  this  pump  is  shown  in  Figure  3.  Pump  design 
ia  baaed  on  a  working  pressure  of  3,  000  psi.  as  a  result,  when  supply  pres¬ 
sures  less  than  3,  000  psi  are  used,  an  over -designed  condition  will  occur  and 
a  weight  penalty  realized.  For  pressures  greater  than  3,  000  psi,  an 
appropriate  correction  factor  is  added  to  the  basic  3,  000  psi  relationship. 


(2)  Accumulator 


Cylindrical  piston  type  and  spherical  bladder  type  accumulators  are  available 
for  this  type  of  application.  Because  Douglas  originally  proposed  the  spheri¬ 
cal  type  for  the  C-5A,  it  will  be  used  in  this  study.  Familiarity,  succeasfui 
utilization  in  past  applications,  and  lighter  weight  are  the  primary  factors 
which  influenced  the  decision  in  favor  of  the  spherical  accumulator.  A 
photograph  of  a  typical  uidt  is  shown  by  Figure  4. 

Because  no  limit  exists  upon  the  duration  of  the  aircraft's  requirement  for 
maximum  power,  the  pump  is  sized  to  accommodate  the  total  required  power 
at  rated  revolutions  per  minute  and  the  accumulator  will  be  used  to  minimize 
the  effect  of  pressure  disturbances  in  the  hydraulic  system. 

(3)  Reservoir 

The  reservoir  will  be  a  boot-strap  type  unit.  It  eliminates  the  adverse  effects 
in  the  hydraulic  system  caused  by  changes  of  aircraft  attitude  and  acceleration 
loadings,  and  because  it  is  airless,  does  not  contribute  to  the  undesirable 
quantity  of  air  dissolved  in  the  hydraulic  fluid.  A  photograph  of  this  type  of 
reservoir  is  shown  by  Figure  5. 

(4)  Tubing 

The  supply  and  return  system  tubing  has  been  divided  into  two  sections. 
Section  No.  1'  (as  shown  in  Figure  2)  is  from  the  pump  to  the  place  where  the 
plumbing  branches  out  to  each  elevator  valve  actuator  package.  Section  No.  2 
includes  all  of  the  tubing  from  the  branch  point  to  Elevators  1  and  2  or  3  and 
4.  Tubing  sizes  can  vary  between  sections  and  between  the  supply  and  return 
systems;  however,  within  a  section,  the  sizes  are  constant.  Steel  tubing  is 
used  for  the  supply  system  and  aluminum  for  the  return  system. 

(5)  Valve  Actuator  Package 

The  valve  actuator  package  is  tailored  after  the  uuit  used  by  Boeing  for  the 
727  elevator -actuation  system.  Figures  6  and  7  are  photographs  of  this  unit. 
The  various  relationships  that  pertain  to  this  unit  are  divided  into  the  follow¬ 
ing  groups:  actuator,  valve  assembly,  and  valve -actuator  package.  Even 
though  the  valve  and  actuator  are  tandem  units,  the  one  hydraulic  system wiil 
be  used  to  supply  both  parts  of  these  units.  Also,  for  simplification,  the 
same  input  signal  will  be  used  for  each  elevator  so  that  the  flow  is  equally 
divided. 

(6)  Mechanical  Linkage 

The  C-5A  flight  control  system  included  a  mechanical  link  between  the  pilot 
and  the  elevator  control  surfaces.  Because  the  weigh*  of  the  mechanical,  link 
was  about  15%  of  all  the  components  and  subassemblies  directly  associated 
with  the  elevator  system,  the  linkage  weight  was  included  in  the  model.  None 
of  the  independent  parameters  influenced  its  design;  therefore,  the  weight 
was  considered  a  constant.  A  schematic  is  shown  by  Figure  8. 
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Figure  7.  Valve  Actuator  Package  -  Too  View 


b.  Dynamic  Model 


The  time -dependent  relationships  which  constitute  the  dynamic  model  are 
concerned  solely  with  the  dynamic  performance  of  the  actuation  system.  The 
development  of  a  function  to  represent  dynamic  performance  is  complicated 
by  the  numerous  quantities  that  must  be  considered  to  adequately  describe  it. 
Most  measures  of  dynamic  performance  are  stated  in  terms  of  the  closed 
loop  behavior  of  the  system  (for  example,  its  bandwidth  and  rise  t'me).  It  is 
particularly  difficult  to  relate  these  terms  to  the  independent  variables  used 
in  the  optimization  procedure.  It  is  clear  that  changing  the  terms  involving 
the  system  pressure,  actuator  area,  and  actuator  moment  arm  in  the  trans¬ 
fer  function  without  changing  the  system  compensation  is  meaningless.  As  a 
result,  it  was  decided  initially  to  use  a  method  to  derive  a  measure  of  dynamic 
performance  which  utilized  control  loop  gains  as  well  as  the  actuation  system 
design  parameters.  The  integrated-square -error  (JSE)  criterion  was  selected 
for  this  purpose  because  a  system  designed  to  this  criterion  results  in  an 
optimal  following  system.  That  is,  the  optimum  system  will  be  defined  as 
the  system  that  yields  a  minimum  integrated-square  error  between  a  refer¬ 
ence  input  and  the  system  output  (see  Figure  9).  Thus,  the  optimum  system 
will  provide  the  best  reproduction  of  the  input  signal.  In  addition,  the  ISE 
was  selected  because  it  provided  a  numeric  that  described  the  optimum  sys¬ 
tem  and  because  a  reference  input  having  the  desired  bandwidth  and  rise  time 
could  be  chosen. 

Because  dynamic  performance  was  being  optimized  for  only  10%  input  signals, 
it  was  possible  to  use  a  linear  representation  of  the  valve -actuator -load  com¬ 
bination  as  shown  by  the  actuation  block  of  Figure  9.  This  linear  model  is 
developed  in  Appendix  II. 

The  ISE  can  be  calculated  analytically  and  the  optimal  parameters  evaluated 
for  systems  of  low  order.  For  this  study,  however,  a  two -level  optimization 
procedure  is  required;  first,  the  system  independent  variables  must  be  speci¬ 
fied  and  then,  second,  the  parameters  of  the  control  system  such  as  Ka  (the 
amplifier  gain)  have  to  be  evaluated.  This  results  in  nonlinear  algebraic 
equations  which  are  hard  to  solve;  therefore,  an  alternate  technique  was 
employed  because  for  the  purpose  of  the  study,  only  the  ISE  is  needed. 

The  alternate  technique  consisted  of  simulating  the  system  on  the  analog  com¬ 
puter  and  evaluating  the  minimum  ISE.  This  was  done  by  first  setting  the 
parameters  of  the  actuation  system  and  then  selecting  the  control  system 
gains  to  yield  the  minimum  ISE. 

From  the  development  of  the  linear  model  included  in  Appendix  II,  it  is  seen 
that  a  change  in  system  pressure  will  affect  only  the  Kq  term.  Because  this 
term  appears  only  in  loop  gains  involving  independent  control  system  para¬ 
meters,  it  is  clear  that  the  system  can  be  made  independent  of  pressure 
changed.  Therefore,  the  oniy  terms  that  can  alfect  the  system  ISE  are  the 
actuator  area  and  the  actuator  moment  arm  thus  making  it  an  implicit  rather 
than  an  explicit  function  of  the  pressure  losses  through  the  plumbing  and  the 
metering  valve.  Further  inspection  of  the  linearized  block  diagram  leads  one 
to  the  conclusion  that  changes  in  area  and  moment  arm  should  be  identical 
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Figure  9.  Integrated  Square  Error  Generation 


within  a  gain  change  in  the  forward  loop.  Simulation  of  the  system  on  the 
computer  verified  this  result  and  allowed  the  ISE  to  be  expressed  as  a  function 
of  the  product  of  area  and  moment  arm.  Because  the  actuator-load  resonant 
frequency  is  also  a  function  of  this  same  product,  it  was  found  that  ISE 
decreased  when  the  product  increased. 

One  disadvantage  of  the  ISE  is  that  it  does  not  indicate  when  the  dynamic  per¬ 
formance  is  adequate.  That  is,  because  a  following  device  is  being  designed, 
any  increase  in  open-loop  resonant  frequency  that  allows  more  faithful  follow 
ing  of  the  reference  input  will  result  in  a  decreased  ISE,  However,  if  the 
previous  system  were  adequate,  the  improved  performance  indicated  by  the 
ISE  is  misleading  because  the  improvement  is  not  required.  Another  disad¬ 
vantage  of  the  ISE  is  that  no  knowledge  is  available  of  how  well  the  system 
could  have  been  made  to  perform  if  additional  and/or  different  types  of 
compensation  had  been  used.  For  these  reasons,  it  was  felt  that  the  ISE 
created  an  artificial  function  and  did  not  contain  adequate  design  information. 

In  an  effort  to  define  a  new  measure  of  dynamic  performance  that  would  be 
more  meaningful  (especially  to  the  designer),  a  function  of  the  open-loop 
resonant  frequency  was  sought. 

Experience  at  Douglas  has  shown  that  it  is  possible  to  achieve  adequate  per¬ 
formance  of  the  control  system  if  the  load-actuator  resonant  frequency  is  at 


least  four  times  the  desired  bandwidth  of  the  actuation  control  loop.  In  some 
instances,  it  was  possible  to  reduce  this  factor  to  as  low  as  2.  5.  Several 
things  should  be  bor*:^  in  mind  concerning  these  statements:  (1)  the  size  of 
this  factor  is  also  a  function  of  the  flight  control  loop  within  wh'ch  the  actua¬ 
tion  system  is  placed;  (2)  this  work  was  done  entirely  with  flow  control  valves 
and  the  load-actuator  poles  are  virtually  undamped  in  this  case;  and  (3)  factors 
greater  than  four  allow  the  designer  greater  flexibility  and  simplify  the  prob¬ 
lem  of  compensation.  Figure  10  shows  a  plot  of  the  open-loop  poles  and  zeros 
of  the  system  transfer  function  as  a  function  of  the  load-actuator  resonant 
frequency.  The  damping  associated  with  these  poles  is  caused  by  the  dynamic 
pressure  feedback.  An  additional  feature  of  a  function  based  only  on  the 
resonant  frequency  is  that  no  compensation  of  the  control  loop  need  be  con¬ 
sidered  as  was  the  case  with  the  ISE.  Referencing  the  block  diagrams  in 
Appendix  II,  the  load -actuation  resonant  frequency  per  actuator  can  be  written 
as 
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The  volume  of  the  actuator  can  be  written  as 
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Fipre  10.  Resonant  Frequency  Root  Locus 


The  resonant  frequency  beumes 


Numerous  function* 
measure  of  dynamic 
Section  IV. 
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SECTION  III 

OPTIMIZATION  TECHNIQUE 


1.  PROBLEM  STATEMENT 

The  problem  under  consideration  is  the  method  for  choosing  opurnal  system 
parameters.  The  optimal  system  is  defined  as  the  system  that  maximizes 
the  chosen  performance  index  (P.  I.  ).  For  this  actuation  design  problem,  the 
system  pressure,  actuator  area,  actuator  moment  arm,  and  tube  diameters 
are  chosen  to  maximize  a  function  of  the  weight,  volume,  dollar  cost,  relia¬ 
bility,  dynamic  performance,  and  environmental  sensitivity.  To  state  this 
mathematic  ally,  let  S  be  a  physical  system  or  process. 


The  parameter  vector  u  is  the  system  pressure,  actuator  area,  actuator 
moment  .arm,  and  tube  diameters  and  x  is  the  output  vector  and  is  composed 
of  the  weight,  volume,  dollar  cost,  reliability,  dynamic  performance,  and 
environmental  sensitivity.  Therefore,  P.  I.  Qp-p(u)  must  be  defined.  Let 

P.I.OPT(u)  =  MAX  P .  I.  (u) 

n 

p.  I.  <u)  =  ^  fjj.(x) 
i  =  1 

Xj  =  g.(u)  j  s  1,  1 . m 

In  general,  the  equations  defining  the  s/stem  will  be  nonlinear.  The  math 
model  derived  in  the  previous  section  related  the  parameter  vector  to  the 
system  outputs.  The  following  sections  relate  the  system  outputs  to  the  cost 
functions,  J.,  and  the  cost  functions  to  the  total  performance  criterion. 

2.  COST  FUNCTION  DEVELOPMENT 

The  development  of  the  cost  functions  is  a  task  that  requires  considerable 
judgment  by  the  designer.  One  fundamental  problem  is  the  generation  of  a 
numeric  that  will,  adequately  represent  the  vehicle  constraints.  Also,  it  is 
desirable  to  have  common  bounds  on  all  of  the  variables.  If  this  is  not  done, 


the  difference  in  the  absolute  values  of  the  different  system  outputs  will  affect 
the  optimization  procedure.  Therefore,  a  numeric  bounded  between  zero  and 
one  was  used  for  the  cost  functions.  Thus,  in  a  maximum  searching  system, 
one  corresponds  to  the  best  possible  system,  and  zero  to  the  worst. 

Develop  ment  of  the  cost  function  to  be  used  on  each  vehicle  constraint  will 
depend  on  the  design  requirements  of  that  vehicle.  For  example,  if  all  incre¬ 
mental  decreases  in  system  weight  are  equally  important,  a  linear  weighting 
on  the  to*al  system  weight  might  be  desirable.  One  method  of  developing  a 
linear  weighting  is  to  compare  the  calculated  weight  with  a  reference  weight. 
For  the  study,  the  reference  weight  is  dei  ned  as  the  weight  calculated  with 
use  ot  standard  design  techniques.  With  use  of  this  cost  function,,  absolute 
changes  from  the  reference  weight  (Wn)  are  weighted  identically.  One 
possible  form  of  the  cost  function  is 


J. 

l 


1 
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n 


Thus,  whe->  the  veight  is  equal  to  the  reference  weight,  J.  «  l/li„  and  varies 
as  shown  by  Figure  11. 

Note  that  the  cost  function  could  be  leso  than  zero  if  the  weight  increased  to 
more  than  twice  the  reference  weight.  If  it  is  mean'ngful  to  consider  such 
cases  for  the  system  being  designed,  the  cost  function  must  be  redefined. 


Figure  11.  Weight  Cost  Function 
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For  other  systems,  any  actuation  system  weight  less  than  a  maximum  allow¬ 
able  weight,  (WmaXi  ),  might  be  permissible.  For  such  a  system,  a  cost 
function  such  as 

might  be  appropriate.  This  cost  function  is  zero  when  W  ,  and  is  one 

when  W  <  Wmax  .  Beside,  the  two  possible  cost  functions  mentioned,  many 
interesting  linear  and  nonlinear  weightings  could  be  derived  to  meet  the 
requirements  of  a  given  design  problem. 

A  vehicle  constraint  that  requires  special  consideration  is  the  volume  of  the 
system.  Consideration  of  the  elements  that  comprise  an  actuation  system, 
leads  to  the  conclusion  that  total  volume  is  not  a  meaningful  quantity  for  many 
of  the  elements.  For  instance,  the  length  of  the  actuator  is  a  more  critical 
quantity  than  the  total  volurr. a  of  the  valve  actuator  assembly  because  the 
valve  can  be  shaped  to  iit  the  space  available.  Therefore,  volume  per  se  was 
not  used  as  the  cost  function.  A  quantity  that  is  defined  in  this  study  as  a 
space  factor  was  substituted.  The  spa..e  factor  is  a  combination  of  each  ele¬ 
ment's  critical  dimensions,  such  as,  actuator  length,  volume  of  the  pump  and 
accumulator,  and  the  length  and  diameter  of  the  reservoir.  One  method  of 
deriving  a  space  factor  is  to  define  ~  subcost  function  for  each  term  of  the 
space  factor.  The  space  factor  is  defined  as  the  sum  of  these  subcost  func¬ 
tions.  With  use  of  the  ideas  presented  earlier,  a  suitable  cost  function  might 
be  defined  as  fallows:  if  the  volume  of  the  elements  for  which  volume  is  the 
important  factor  is  summed,  a  suitable  subgoal  might  be 

V,  =  1  -  Veq/2  Veq 
1  mom 

where  Veq  is  the  current  volume  and  Veqnom  is  the  refc :ence  volume  which 
might  be  the  value  given  by  standard  design  techniques.  If  the  length  of  the 
actuator  is  a  critical  dimension,  a  suitable  subgoal  might  be 

V,  =  1  -  ACTL/ZACTL 
2  nom 

where  ACTL  is  the  current  actuator  length  and  ACTLnom  is  a  reference 
length.  The  space  factor  cost  function  might  be  defined  as 

J  =  1  -  VOL/2VOL 

nom 

where  VOL  is  the  current  value  of 

n 

VOL  =  - - - Y  f.V. 

r-  y  f.  ^ 
i  =  1 
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and  VOLjj0m  is  a  value  calculated  as  a  reference.  The  f.'s  are  emphasis 
factors  which  must  be  specified  by  the  designer. 

Each  critical  dimension  should  be  assigned  a  subgoal  and  the  weighting 
coefficients  chosen  to  represent  their  importance  to  the  total  cost  function. 

The  dollar  cost  function  exhibited  no  unusual  characteristic  and  was  carried 
along  as  a  linear  cost  function  in  this  study. 

The  reliability  cost  function  was  found  to  be  almost  independent  of  the  inde¬ 
pendent  variables  of  the  study  and  was  simply  carried  along  as  a  linear  cost 
function. 

The  environmental  sensitivity  cost  function  is  another  factor  that  requires 
some  special  attention.  It  is  clear  that  dynamic  performance  and  environ¬ 
mental  sensitivity  are  closely  related.  In  terms  of  the  actuator  design  problem, 
the  only  term  in  the  environmental  sensitivity  cost  function  would  be  the  bulk 
modulus  of  the  oil.  This  term  directly  affects  the  actuator  load  resonant  fre¬ 
quency,  which  is  the  term  that  affects  the  dynamic  performance  cost  function. 
A  definition  of  an  environmental  sensitivity  cost  function  that  might  be  mean¬ 
ingful  is  the  comparison  of  the  worst-case  dynamic  performance  with  the 
nominal  dynamic  performance  for  each  design.  Thu3,  if  max.  (Dyn.  P.  )  is 
defined  to  be  the  maximum  (or  minimum)  value  of  the  dynamic  performance 
cost  function  caused  by  the  change  in  bulk  modulus  of  the  oil,  then 


max.  (Dyn.  P.  )  -  Dyn.  P. 
- 


This  function  could  exceed  the  bounds  of  one  if  the  maximum  value  of  dynamic 
performance  were  greater  than  twice  the  nominal  value.  This  would  indicate 
extreme  sensitivity  to  environmental  factors  and  that  the  design  should 
probably  be  rejected. 

3.  SYSTEM  PERFORMANCE  INDEX 

As  noted  earlier,  a  performance  index  that  can  be  us  sd  to  indicate  the  total 
merit  of  a  particular  design  is: 


n 


i  «1 


where  the  J^'s  again  are  the  cost  functions.  The  f^  terms  are  called  emphasis 
coefficients  and  are  used  to  denote  the  particular  importance  ot  each  cost 
function  to  the  design.  If  '.e  make  the  followi~j  restriction 


Ih 


100 
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then  the  fj’s  can  be  used  to  denote  a  percentage  importance  of  each  cost  func¬ 
tion  to  the  whole  system.  It  is  apparent  why  the  J;'s  were  bounded  between 
zero  and  one;  otherwise  the  emphasis  coefficients  would  lose  their  connota¬ 
tion  of  percentage  emphasis. 

4.  COMPUTER  PROGRAMS 

The  computer  programs  developed  for  the  optimal  design  problem  are  tool? 
for  the  evaluation  of  the  performance  indexes.  Two  separate  programs  were 
developed  to  provide  information  to  the  designer. 

The  first  program  (No.  1)  was  developed  primarily  f_r  providing  information 
concerning  the  sensitivity  of  the  cost  functions.  This  is  important  for  choos¬ 
ing  the  emphasis  coefficients  because  the  a  priori  guess  may  lead  to  a  domi¬ 
nation  of  the  P.  I.  by  one  cost  function.  TRe  maximum  deviation  of  each  cost 
function  must  be  considered  in  context  v/ith  its  emphasis  coefficient.  To  pro¬ 
vide  the  sensitivity  information,  the  parameter  space  was  systematically 
searched  by  imposing  a  fixed  grid  on  the  space.  The  cost  functions  and  the 
performance  index  were  evaluated  and  printed  at  each  point  on  the  grid.  This 
information  can  be  evaluated  and  the  emphasis  coefficients  updated  as  required. 
A  -low  diagram  representing  the  essential  features  of  this  program  is  shown 
in  Figure  12.  A  detailed  description  and  a  listing  of  the  program  are  located 
in  Appendix  III. 

The  second  program  (No.  2)  developed  was  primarily  intended  for  choosing 
the  optimal  system.  Nun.erous  -echniques  of  linear  and  nonlinear  program¬ 
ming  are  available  for  this  purpose.  Because  there  are  bounds  on  the  param¬ 
eter  space  and  on  the  cost  functions,  the  problem  is  fundamentally  a  nonlinea." 
programming  problem.  One  of  the  more  popular  searching  techniques  is  the 
gradient  method.  Application  of  this  technique  to  this  problem,  however,  is 
complicated  by  the  inclusion  of  the  tube  diameters  as  independent  variables. 
The  tube  diameters  are  discrete  variables;  that  is,  they  take  on  only  a  finite 
number  of  states,  and  special  techniques  must  be  used  for  handling  discrete 
variables.  Another  popular  searching  technique  is  the  random  search  or  an 
extension,  the  adaptive  random  search.  Its  chief  advantage  is  the  simplicity 
of  the  programming.  It  has  the  disadvantage,  however,  that  it  may  take 
excessive  run  time  on  the  computer.  There  are  two  tradeoffs  to  consider: 

(1)  whether  it  is  more  efficient,  in  terms  of  the  computer,  to  calculate  the 
gradient  or  to  accept  the  wrong  guesses  of  the  random  search,  or  (2)  whether 
it  is  more  efficient,  in  terms  of  manpower,  to  use  a  simple  program  at  the 
expense  of  computer  time.  After  some  experience  with  the  fixed-grid 
searching  program,  it  was  found  that  the  program  took  little  machine  time, 
so  the  adaptive  random  search  program  was  chosen.  A  flow  diagram  repre¬ 
senting  the  essential  features  of  this  program  is  shown  in  Figure  13.  Detailed 
descriptions  of  the  random  search  technique  and  of  the  No.  2  program  are 
located  in  Appendix  IV.  Listings  that  differ  from  the  No.  1  program  are  also 
included. 
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SECTION  IV 
APPLICATION 


The  cost  functions  for  the  sample  problem  are  developed  in  this  section  and 
the  results  of  the  optimization  procedure  are  discussed.  The  design  require¬ 
ments  and  restrictions  are  given  in  paragraphs  2  and  3  of  Section  II.  To 
supply  a  reference  value  for  all  cost  functions  using  a  linear  weighting,  a 
nominal  system  was  designed  using  the  design  criterion  of  the  Douglas  pro¬ 
posal  for  the  C-5A.  This  also  provided  a  convenient  measure  of  the  merit  of 
the  new  design. 

1.  COST  FUNCTIONS 

a.  Weight  Cost  Function 

In  the  derivation  of  the  sample  problem,  it  was  felt  that  weight  in  each  por¬ 
tion  of  the  actuation  system  was  equally  important;  in  addition,  it  was  felt 
that  increases  and  decreases  in  the  weight  of  each  portion  of  the  system  were 
also  equally  important.  Therefore,  a  cost  function  using  the  total  weight  of 
the  system  with  a  linear  weighting  was  used.  That  is 

J.  =  1  -  W/2W 
l  norn 

When  it  was  felt  that  weight  in  one  section  of  the  system  was  more  important 
than  weight  in  another  section,  a  compound  cost  function,  similar  to  the 
volume  cost  function  described  earlier,  could  be  used.  F->r  example,  weight 
that  is  located  a  considerable  distance  from  the  center-of-gravity  (CG)  of  the 
aircraft  might  be  more  critical  than  weight  located  near  the  CG  and,  hence, 
the  former  weight  should  receive  more  emphasis  in  the  cost  function  than  the 
latter . 

b.  Space  Factor  Cost  Function 

With  use  of  the  ideas  presented  in  Section  III,  the  space  factor  cost  function 
was  defined  as  a  sum  of  the  critical  dimensions  of  the  elements  of  the  actua¬ 
tion  system.  An  equivalent  volume  was  defined  as  the  sum  of  the  volumes  of 
the  elements  for  which  total  volume  was  thought  to  be  an  important  factor. 

For  this  example,  only  the  pump  and  the  accumulator  were  included  in  this 
sum.  Equivalent  volume  was  added  as  one  of  the  subgoals  of  the  space  factor 
cost  function.  The  reservoir  was  not  a  shape  that  could  be  adequately 
described  by  a  volume.  It  was  thought  that  the  diameter  of  the  lc-ge  end  and 
the  length  of  the  reservoir  were  a  better  description.  Therefore,  these  two 
critical  dimensions  were  each  included  as  subgoals.  Another  critical  dimen¬ 
sion  was  the  length  of  the  actuator.  It  was  felt  that  volume  was  not  an  ade¬ 
quate  descriptor  because  the  valve  can  be  shaped  to  fit  the  volume  available. 
Length,  howevei ,  can  be  very  critical  in  terms  of  the  design  of  the  aircraft 
structure.  For  this  problem,  the  diameter  of  the  actuator  was  not  critical. 
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In  any  aircraft  with  a  narrow  section,  however,  it  might  be  very  important. 
The  original  definition  cf  the  cost  function  included  these  four  terms  with 
equal  weighting  and  with  a  linear  subgoal  for  each  term.  The  reference  value 
for  each  of  the  linear  subgoals  was  the  value  obtained  using  standard  design 
practice.  Thus,  the  cost  function  is  written  as 


*  -r-  t  ftQi 


i“  1 


=  1 


Where 


Q,  =  1  - 


Veq 


10m 


Q  =  I  .  — actl 
U2  2ACTL, 


nom 


Q  -  i  .—REP 

U3  2RED 


nom 


o  -  l  _ . . .  .RBilrf-i 

U4  "  1  2REL 


nom 


In  the  design  of  the  C-5A  structure,  the  space  available  for  the  actuator  was 
limited  in  length  by  the  supporting  structure.  To  include  the  effect  of  this 
design  restriction,  an  additional  weighting  was  placed  on  any  design  that 
required  a  length  of  actuator  longer  than  the  length  readily  available.  That 
is,  a  term 


1/8  (1  -  sgn  (5.5-R)} 


was  added  to  the  actuator  length  subgoal,  where  5.5  in.  is  the  maximum 
length  readily  available.  The  term  penalizes  the  use  of  any  actuator  longer, 
than  the  maximum  length  on  an  equal  basis.  This  restriction  is  not  meant 
to  limit  the  length,  but  to  give  some  importance  to  the  extra  problems 
involved  in  using  the  longer  moment  arm. 

n  -ff  Co?t  Functlor. 

As  discussed  in  Appendix  II,  two  variations  were  investigated.  The  ISE 
model  indicates  improved  performance  whenever  the  resonant  frequency  of 
the  load-actuator  combination  is  increased.  This  presents  a  problem, 
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because  system  requirements  already  have  been  established.  Investigation 
of  the  system  on  the  analog  computer  showed  that  performance  requirements 
of  the  system  were  met  easily  because  of  the  high  hinge-moment  require¬ 
ments .  It  w?s,  therefore,  felt  that  the  ISE  was  misleading,  because  the  sys¬ 
tem  was  not  as  critical  as  the  ISE  indicated.  For  this  reason,  plus  the  fact 
that  designers  are  not  normally  familiar  with  ISE,  an  alternate  definition  of 
dynamic  performance  was  selected.  The  alternate  selection,  which  is  in 
terms  of  system  bandwidth,  is 

TOTDYP  =  1  -  EXP  [-l/4(u  -  O.lu)] 

and  is  shown  in  Figure  14.  TOTDYP  is  the  output  of  the  d>namic  perfor¬ 
mance  subroutine.  The  cost  function  is  one  minus  TOTDYP. 

There  ts  a  good  reason  for  choosing  this  cost  function.  As  stated  in  the  math 
model  of  the  system,  experience  shows  that  it  is  possible  to  compensate  the 
system  if  the  open  loop  poles  are  at  least  four  times  farther  out  in  frequency 
than  the  desired  bandwidth.  The  cost  function  places  no  penalty  on  the  sys¬ 
tem  if  the  resonant  frequency  is  10  times  farther  oul  in  frequency  than  the 
desired  bandwidth.  As  the  resonant  frequency  decreases,  an  increasing 
penalty  is  placed  on  the  design.  Because  the  cost  function  could  become 
negative  for  frequencies  greater  than  10  times  the  bandwidth,  the  function 


Figure  14.  Dynamic  Performance  Function 
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was  defined  to  be  zero  for  any  value  of  the  exponential  '  -enter  than  one. 

Thin  indicates  that  there  is  no  particular  value  in  havin'  a  wider  bandwidth 
than  required. 

The  cost  function  shown  in  Figure  14  may  not  be  desirable  because  it  heavily 
weights  small  decreases  in  resonant  frequency.  A  more  suitable  cost  func¬ 
tion  is  shown  in  Figure  15. 

This  cost  function  weights  small  changes  very  lightly,  but  the  weighting  gets 
prog  ressively  heavier  until  it  reaches  the  maximum  at  4uq. 

It  is  possible,  in  fact,  that  this  excess  bandwidth  could  be  a  problem  in  terms 
of  bending  modes  or  flutter  frequencies  of  the  aircraft.  If  certain  frequen¬ 
cies  must  be  avoided,  it  can  be  handled  easily  with  this  type  of  cost  function 
by  adding  another  term  that  takes  a  value  of  one  at  the  frequencies  in  ques¬ 
tion,  Such  a  term  might  be 


EXP  (uq  -  u)2J 

This  term  has  the  desired  effect  of  heavily  weighting  frequencies  close  to  the 
critical  frequency  and  not  weighting  the  other  frequencies  at  ail. 


FREQUENCY 


Figure  15.  Alternate  Dynamic  Performance  Function 


One  additional  point  should  be  noted;  an  additional  refinement  of  dynamic 
performance  could  be  made  by  including  the  functional  relationship  between 
system  elasticity  and  the  system  pressure  and  area.  For  the  example  in  this 
report,  the  effect  was  not  significant  because  the  resonant  frequency  was 
very  high  and  the  change  in  system  elasticity  is  not  great  at  the  temperature 
and  pressures  involved.  In  problems  where  the  resonant  frequency  is  criti¬ 
cal,  this  maybe  an  important  factor. 

d.  Environmental  Sensitivity  Cost  Function 

3ecause  the  elasticity  of  the  system  directly  influences  the  resonant  fre¬ 
quency  of  the  system,  the  environmental  sensitivity  cost  function  is  defined 
easily  in  terms  of  the.  dynamic-performance  cost  function.  For  the  example 
problem,  the  bulk  modulus  of  the  oil  and  the  elasticity  of  the  actuator  struc¬ 
ture  combine  to  change  the  elasticity  of  the  system.  However,  it  was  found 
that  this  effect  was  not  significant,  especially  because  the  dynamic  perfor¬ 
mance  was  not  a  significant  factor.  For  this  reason,  this  co*t  function  was 
not  considered  further . 

e.  Dollar  Cost  Function 

The  dollar  cost  function  was  taken  to  be  an  increment?!  cost  and  was  not 
strongly  a  function  of  the  independent  variables.  Because  of  its  lack  of 
significance,  it  was  taken  to  be  a  linear-cost  function  and  was  assigned  a 
small  weighting. 

f.  Reliability  Cost  Function 

Although  it  was  felt  that  reliability  would  play  a  key  role  in  the  choice  of  an 
optimum  system,  sufficient  data  could  not  be  found  to  establish  a  connection 
between  reliability  and  the  independent  parame‘ers.  The  cost  function  was 
carried  along  as  a  linear  function  with  small  r  phasis,  in  the  event  that 
such  a  significance  could  be  found. 

2.  EMPHASIS  COEFFICIENT^ 


The  relative  importance  of  the  various  vehicle  constraints  is  incorporated 
into  the  optimum  design  by  selection  of  the  emphasis  coefficients.  In  this 
particular  problem,  weight  was  considered  to  be  the  most  important  factor. 
Dynamic  performance  and  volume  were  given  equal  ratings  followed  by  cost 
and  then  reliability.  Not  only  should  the  relative  importance  be  considered 
when  selecting  these  coefficients  but  the  extent  to  which  the  corresponding 
analytical  relationships  are  developed  and  the  accuracy  or  confidence  in  the 
basic  input  data  should  also  be  taken  into  consideration.  The  percentage 
values  selected  are  as  follows: 


1 . 
2. 
3. 


f  .  ,  =  30%. 

weight 

^dynamic  performance 


volume  =  25%. 


25%. 


reliability 

3.  OPTIMIZATION  PROCESS 

With  the  cost  functions  defined  and  the  model  defined  as  in  Section  2,  tha 
problem  was  put  on  the  digital  computer  for  evaluation  of  the  optimum  sys¬ 
tem.  The  computer  program  using  the  fixed-grid  searching  technique  was 
used  to  generate  data  that  could  be  used  by  the  designer  in  evaluating  the 
results.  The  designer  has  to  know  which  cost  functions  are  the  most  depen¬ 
dent  on  the  independent  variables.  This  is  important  for  the  proper  choice  of 
the  emphasis  coefficients  and  for  an  understanding  of  the  results  of  the  opti¬ 
mization  process.  The  program  pointed  up  the  problem  in  the  use  of  ISE.  It 
was  found  tba*  only  the  ISE  showed  a  significant  change  with  the  independent 
variables  and  that  the  dynamic-performance  cost  function  dominated  the  per¬ 
formance  index.  A  review  of  this  problem  led  to  a  redefinition  of  the 
dynamic -performance  cost  function  which  is  felt  to  more  adequately  describe 
the  physical  situation. 

All  of  the  other  cost  functions  showed  little  variation  with  changes  in  the 
independent  variables.  As  can  be  seen  in  Figures  16  through  21,  the  graphs 
have  large  nominal  values  and  show  little  curvature.  Figure  1 6  shows  a 
gra  <h  of  weight  as  a  function  of  pressure  for  a  constant-torque  system;  that 
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Figure  Ib.Weight  as  a  Function  of  Pressure 
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Figure  21.  P.l.  as  a  Function  of  My...*nt  Arm 


is,  at  each  pressure,  the  minimum  actuator  area  needed  to  supply  the  hinge 
moment  and  control  surface  rate  for  the  given  moment  arm  is  used.  The 
minimum  area  is  determined  by  calculating  that  flow  rate  and  corresponding 
actuator  AP  necessary  to  meet  requirements  taking  into  consideration  the 
pressure  losses  through  the  plumbm-'  and  the  metering  valve.  This  curve 
shows  a  slight  decrease  as  pressure  increases  from  2,000  to  3,000  psi. 

This  decrease  is  primarily  because  cf  the  smaller  items  required  at  the 
higher  pressures.  Above  3,  000  psi,  this  decrease  is  lost  in  the  inc:  sase  in 
weight  caused  by  the  weight  penalty  to  maintain  the  required  safety  factor  . 
The  minimum  weight  occurs  at  3,  000  psi  and  is  745  lb.  Note  that  this  mini¬ 
mum  at  3,  000  psi  is  somewhat  artificial  because  of  the  weight  penalty  abovt 
3,  000  psi  and  the  fact  that  below  3,  000  psi,  standard  3,  000  psi  equipment  is 
used,  If  the  2,000  ps;  equipment  were  redesigned  to  reflect  the  lowered 
stresses  there  might  not  be  a  minimum  at  3,  000  psi.  This  is  not  standard 
practice,  however.  Disregarding  this  fact,  the  curve  chows  that  most  of  the 
actuation  system  weight  is  fixed  and  that  the  pressure  can  be  chosen  over  a 
wide  ’•ange  with  little  penalty  in  weight.  For  example,  the  maximum  devia- 
t’on  in  weight  between  2.  000  and  6.  000  psi  is  approximately  9. 3%.  or 
appi  oximately  70  lb.  It  is  probable  that  such  a  small  difference  would  not 
be  significant  in  a  large  aircraft  system. 

Figure  17  shows  a  graph  of  sptce  factor  as  a  function  of  pressure  for  con¬ 
stant  torque  ana  for  a  moment  arm  of  5.5  in.  The  graph  shows  that  the 
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epace  factor  (before  normalization)  decreases  monotonically  from  a  high  of 
2.20  at  2,  000  psi  to  a  low  of  1.69  at  6,  000  psi .  The  r  pace  factor  decreases 
because  the  equivalent  volume  is  a  function  of  flow  which  decreases  as  the 
pressure  increases. 

The  dollar  cost  of  the  system  varied  from  a  high  of  $10,  250  at  2,  000  psi 
monotonically  to  a  low  of  $10,  034  at  a  pressure  of  6,  000  psi.  The  difference 
in  these  numbers  is  lees  than  the  confidence  present  in  the  values  them¬ 
selves.  The  maximum  deviation  is  approximately  0.2%. 

Reliability  was  found  to  be  independent  of  pressure. 

When  the  revised  dynamic  performance  index  was  used,  dynamic  perfor¬ 
mance  was  found  to  be  independent  of  pressure. 

Figure  18  shows  a  graph  of  weight  as  a  function  of  the  moment  arm  at  a 
pressure  of  3,  000  psi  and  at  constant  torque.  This  shows  that  weight  is 
virtually  indepsndent  of  the  moment  arm--a  maximum  deviation  of  less  than 
0.2%.  The  graph  of  actuator  weight  as  a  function  of  the  moment  arm  in 
Appendix  I  (Figure  41)  shows  marked  curvature,  but  the  difference  in  scales 
is  important.  One  point  that  might  influence  the  weight  as  a  function  of 
moment  arm  that  was  not  included,  is  the  increased  weight  of  the  supporting 
structure  as  the  arm  decreases.  The  support  was  considered  a  rigid  body  in 
the  study. 

The  space  factor  as  a  function  of  moment  arm  ot  a  constant  torque  for  a 
pressure  of  3,000  psi  is  shown  in  Fig  ire  19.  The  discontinuity  reflects  the 
additional  weighting  for  moment  arms  that  exceed  the  length  that  will  easily 
fit  in  the  available  space.  Again,  the  space  factor  has  not  yet  been 
normalized. 

The  graphs  of  dollar  cost  and  reliability  as  a  function  of  moment  arm  are  not 
shown  because  they  are  not  significant.  Dynamic  performance  is  not  shown 
bee:  use  it  is  always  acceptable. 

The  curves  showing  independent  increases  in  the  actuator  area  aro  not  shown 
because  all  of  the  cost  functions  showed  poorer  performance.  These  data 
were  taken  at  fixed  pressure  and  moment  arm  and,  thus,  reflected  a  system 
of  higher  hinge-moment  capability. 

Figures  20  and  21  show  graphs  of  performance  index  as  a  function  of  pres¬ 
sure  and  moment  arm.  The  P.I.  as  a  function  of  pressure  peaks  at 
4,  000  psi.  The  P.I.  as  a  function  of  moment  arm  is  dominated  by  the  space 
factor  cost  function  and  indicates  the  choice  of  a  short  moment  arm.  The 
effect  of  changing  the  tubing  sizes  is  to  vertically  shift  the  curves;  the 
smaller  the  tube  size,  the  larger  the  performance  index. 

Trie  maximum  f.  1.  is  found  at  4,000  psi,  a  4-j.n.  moment  arm,  minimum 
area,  and  all  of  the  smallest  possible  tuDe  sizes.  At  the  maximum,  the 
value  of  the  P.I.  was  67.49,  the  total  weight  was  638.9  lb,  and  the  space 
factor  was  1.57.  Comparison  with  the  nominal  design  (constant  torque, 
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3,000  psi,  5,5  in.  moment  arm  and  largest  possible  tube  sizes)  indicates  a 
weight  sayings  of  approximately  105  lb  and  a  decrease  in  space  factor  of  0  43, 
An  additional  facet  of  decreased  tube  sizes  that  should  be  noted  is  that  many 
of  the  low  pressures  are  eliminated  by  the  flow  requirement  and  plumbing 
pressure  losses  at  -40°  F. 

Referenced  to  the  maximum  P.I.  using  all  the  larger  tube  sizes  (1  in.  O.D. 
in  Section  1  and  3/4  in.  O.  D.  in  Section  2)  there  is  a  possible  weight  savings 
of  approximately  130  lb  and  a  decrease  in  the  space  factor  cf  0.25  by  uoing 
the  smaller  tube  sizes  (3/4  in.  O.D.  in  Section  1  and  5/8  in.  O.D.  in  Sec¬ 
tion  2).  At.  these  tube  sizes,  the  plumbing  pressure  loss  was  740  psi  for  a 
rated  flow  of  26  gpm  at  a  fluid  temperature  of  120®F.  Use  of  even  smaller 
tube  sizes  was  not  possible  because  of  the  -40°  F  flow  requirement.  The 
above  savings  can  be  realized  if  proper  valve  operation  is  possible  at  lower 
han  normal  valve  pressure  drops. 

Most  of  the  results  presented  are  reasonable,  from  the  point  of  design  prac¬ 
tice.  One  possible  disturbing  item  is  the  increase  in  P.I.  as  the  moment 
arm  decreases.  Possibly,  if  the  weight  of  the  supporting  structure  had  been 
included,  this  effect  would  not  have  been  seen.  Also,  a  short  moment  arm 
increases  the  effect  of  valve-actuator  nonlinearities  on  control  surface 
motion.  To  include  the  latter  it  would  be  necessary  to  make  dynamic  perfor¬ 
mance  a  compound  cost  function,  like  the  space  factor. 

In  an  :ffort  to  gain  additional  insight  into  the  nature  of  the  function  being 
maximized,  data  were  taken  along  constant  pressure,  area,  and  moment  arm 
planes.  Also,  these  data  were  taken  with  the  ISE  dynamic  performance  cost 
function  to  show  what  the  curves  might  be  if  dynamic  performance  was  a 
critical  factor. 

Figure  22  shows  the  graph  of  weight  as  a  function  of  pressure.  It  was  found 
that  space  factor,  ISE,  and  cost  were  constant  as  a  function  of  pressure, 
when  area  and  moment  arm  were  held  constant. 

Figure  23  shows  how  weight,  ISE,  space  factor,  and  coat  vary  as  a  function 
of  area  for  fixed  pressure  and  moment  arm.  The  abscissa  should  be  inter¬ 
preted  as  follows:  (1)  indicates  the  minimum  area  for  the  pressure  and 
moment  arm  given,  (2),  (3),  and  (4)  indicate  increases  of  0.5  in.  ^  of  area 
to  the  minimum  for  each  number.  Note  that  ISE  improves  as  all  of  the  other 
cost  functions  decrease.  Figure  24  shows  the  same  information  for  changes 
in  moment  arm. 

For  comparison,  the  previous  weight  and  cost  data  taken  on  a  cons. ant- 
torque  basis  and  plotted  as  a  function  of  pressure  will  be  shown  to  the  same 
scale  as  used  by  Figures  22  and  23.  Figures  25  and  26  show  these  graphs. 
Note  that  the  scales  on  the  ordinate  are  broken  to  shew  the  curvature.  The 
reason  the  weight  is  not  identical  to  Figure  16  is  becaust  the  linkage  weight 
(which  is  constant)  had  not  yet  been  added  to  the  program  and  only  one  valve 
actuator  assembly  was  being  considered. 

A  similar  comparison  of  cost  functions  as  a  function  of  area  is  not  possible 
because  an  increased  area  reflects  an  increased  torque  capability.  A 
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comparison  of  graphs,  showing  changes  in  the  cost  functions  as  a  function  of 
moment  arm,  is  not  shown  because  the  constant-torque  curves  contained  no 
significant  information. 

The  results  of  this  work  showed  that  tttere  was  little  to  be  gained  in  terms  of 
total  performance  by  the  choice  of  the  independent  parameters  in  this  study. 

4.  RESULTS  OF  THE  RANDOM-SEARCH  TECHNIQUE 

After  the  results  of  the  fixed-grid  search  have  been  evaluated,  it  is  possible 
to  use  the  ranootn- search  routine.  The  fixed-grid  routine  provides  informa¬ 
tion  ab out  the  Bersitivity  of  the  cost  functions  necessary  for  the  choice  of 
emphasis  co*ffic iente  Once  the  emphasic  coefficients  have  been  chosen, 
application  of  the  random-search  routine  is  strsigntforward.  For  this  par¬ 
ticular  problem,  the  random-search  routine  provided  little  information  that 
not  readily  *v»ilahl»  bv  locking  at  results  of  the  fixed-grid  routine.  In 
some  cases,  however,  the  cost  functions  may  vary  a  great  deal  us  a  function 
of  the  independent  variables,  and  the  cost  functions  rnay  be  difficult  to  inter¬ 
pret  from  the  fixed-grid  data.  In  these  cases,  the  optimum  system  will  be 
chosen  by  the  random-searen  procedure.  The  search  procedure  could  also 
give  information  about  the  sensitivity  cf  the  optimum,  if  ail  of  the  successful 
steps  were  printed.  For  the  sample  problem,  it  was  found  that  the 
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random-search  procedure  took  leas  time  to  reach  the  optimum  thaa  the 
fixed-grid  routine.  The  search  procedure  was  terminated  when  the  computer 
made  1,000  unsuccessful  steps  in  succession.  One  problem  that  led  to  the 
high  number  of  unsuccessful  trials  was  the  choice  of  tube  sizes.  The  perfor¬ 
mance  index  was  maximized  when  all  of  the  tube  sizeo  were  at  their  smallest 
value.  Thus,  the  program  had  to  choose  all  of  the  small  tube  sizes,  which  is 
a  low- probability  event  (1/16),  and  also  choose  a  set  of  parameters  (pres¬ 
sure,  area,  and  moment  arm)  that  led  to  a  larger  performance  index.  The 
average  run  took  lese  than  3  min.  of  computer  time  on  the  7094.  This  time 
included  approximately  1  min  of  compilation  time. 

The  adaptive  random- search  routine  showed  little  advantage  for  this  particu¬ 
lar  example.  The  problem  just  ciied,  concerning  the  tube  sizes,  caused  too 
many  unsuccessful  steps  and  negated  the  accumulation  of  knowledge  gained  by 
the  successful  steps.  The  adaptive  search  routine  hae  been  shown  to  be 
successful  in  reducing  search  time  in  other  problems  studied  at  Douglas 
(References  1  and  2). 


SECTION  V 


EVALUATION 


1.  GENERAL 

The  main  objectives  of  this  study  were  to  develop  ?n  optimal  technique  for 
the  design  of  an  actuation  system,  to  determine  guidelines  for  its  use,  and  to 
evaluate  the  usefulness  of  the  technique  with  respect  to  current  design 
methods.  The  results  indicate  that  an  optimal  technique,  and  its  corre¬ 
sponding  guidelines,  have  been  developed.  It  appears,  however,  that  the  data 
obtained  by  application  to  the  sample  problem  are  not  adequate  to  evaluate  its 
merit.  Detailed  discussion  of  the  optimal  technique  and  corresponding  guida- 
lines  is  included  in  Sections  III  and  IV.  A  detailed  evaluation  of  its  useful¬ 
ness  follows. 

2.  TECHNIQUE  EVALUATION 

Whether  evaluating  an  optimal  design  technique  or  any  new  method,  the 
questions  that  must  be  answered  include  the  following: 

a.  Are  the  results  better  than  those  obtained  by  use  of  previous  meth¬ 
ods  and,  if  so,  can  definite  conclusions  be  made  concerning  the 
merit  of  the  new  method? 

b.  How  universal  is  the  new  method,  and  are  the  results  obtained  by 
application  to  different  systems  directly  comparable? 

c.  How  does  the  cost  in  time  and  manpower  compare  between  the  old 
and  new  methods? 

d.  Are  there  related  efforts  being  conducted  elsewhere? 

An  evaluation  by  the  review'  council  provided  answers  that  were  most  appli¬ 
cable  to  the  first  question  because  their  primary  function  was  a  comparison 
between  the  optimal  technique  and  previous  design  methods.  In  general,  this 
evaluation  emphasized  the  fact  that  the  optimal  design  of  the  sample  problem 
did  not  show  much  improvement  o'  er  that  obtained  by  using  current  design 
methods  and  as  a  result,  .t  is  not  possible  to  draw  definite  conclusions  on  the 
merit  of  the  optimal  technique.  The  sample  problem,  a  relatively  elemen¬ 
tary  hydraulic  system,  was  defined  so  as  to  expedite  development  of  the 
technique,  and  so  this  is  possibly  a  reason  for  the  uneventful  results.  It  was 
not  apparent  at  the  beginning  of  the  study  that  more  meaningful  results  would 
not  he  obtained.  As  it  turned  out,  however,  the  relationships  between 
selected  design  parameters  and  vehicle  constraints  were  weak;  so  deviations 
from  nominal  values  were  relatively  insignificant.  If  the  restrictions  and 
boundary  conditinna  the  sample  problem  were  relaxed  and  a  Kieaicr  nu.u- 
ber  of  interrelated  parameters  were  considered,  possibly  more  meaningful 
results  would  be  obtained.  Therefore,  it  fol  -ws  logically  that  this  technique 
should  be  studied  further  by  conducting  addit  nal  tests  utilizing  a  complete 
actuation  system. 


45 


Use  of  an  optimal  technique  with  a  computer  does  not  necessarily  mean  that 
an  answer  will  be  obtained  that  could  not  be  obtained  by  using  current  design 
methods.  On  the  other  hand,  it  does  mean  that  the  monotonous  repetition  of 
conducting  tradeoff  studies  between  the  various  designs  and  design  concepts 
can  be  largely  eliminated.  Also,  by  use  of  this  new  technique,  more  situa¬ 
tions  and  parameter  variations  can  be  investigated,  with  the  hopeful  result 
that  there  is  a  more  expedient  method  to  find  the  optimal  design.  Even 
though  this  optimal  design  is  selected  automatically  by  application  of  the  sys¬ 
tematic  searching  routine  (described  previously) ,  the  results  are  still  subject 
to  the  designer's  background  and  experience  through  the  mathematical  model 
and  other  computer  inputs.  To  standardi-e  the  optimal  design  technique,  it 
is  essential  that  model  assumptions  and  constraint  emphasis  be  controlled  by 
the  agency  conducting  the  design  investigation. 

A  possible  use  of  this  design  technique  is  as  a  synthesis  tool  in  an  overall 
computer  simulation  of  a  complete  vehicle  flight-control  system.  By  doing 
so,  detailed  actuation  system  design  parameters,  not  normally  included  in 
the  initial  design  studies,  could  be  considered. 

The  second  question  is  concerned  with  performance  index  comparison  when 
optimizing  designs  have  different  basic  concepts.  Also,  how  universal  is  the 
technique  itself  when  optimizing  the  eame  systems  which  have  varying 
amounts  of  complexity  and  when  optimizing  systems  which  have  differing 
design  concepts?  In  other  words,  how  much  work  is  involved  when  convert¬ 
ing  this  design  method  for  use  on  different  systems?  Many  such  comparisons 
and  conversions  would  occur  when  conducting  tradeoff  studies.  Because  only 
one  system  was  investigated,  additional  studies  will  be  required  before  these 
questions  can  be  answered. 

When  an  optimal  design  is  determined,  using  this  technique,  the  third  question 
is:  How  much  effort  is  involved  and  how  does  it  compare  with  efforts  required 
in  current  design  methods?  In  an  aerospace  vehicle,  where  many  design 
parameters  are  considered,  use  of  the  synthesis  tool  described  in  this  report 
may  well  be  worth  the  effort  expended.  These  types  of  systems  will  have  to 
be  investigated  further  before  an  answer  can  be  determined. 

Actuation  systems  from  other  types  of  aerospace  vehicles  should  be  included 
in  any  such  investigations.  In  high-performance  missiL*  systems,  for 
example,  dynamic  response  and  power  generation  tradeoffs  have  a  different 
order  of  importance  and  will  have  a  different  influence  on  the  optimal  design. 
As  far  as  the  optimization  procedure  is  concerned,  the  primary  effect  of 
considering  different  aerospace  vehicles  will  be  in  the  choice  of  emphasis 
coefficients.  Cost-function  revisions  should  be  secondary  because  of  in 
attempt  to  build  this  generality  into  the  basic  program.  Of  course,  a  new 
set  of  equations  relating  the  parameter  vector  to  the  output  vector  must  be 
derived. 

The  fourth  question  is  concerned  with  related  efforts  in  the  aerospace  indus¬ 
try,  of  which  the  principal  one  is  the  study  being  conducted  by  General 
Electric,  which  has  an  objective  similar  to  this  study’s.  However,  the  work 
reported  in  this  report  deals  with  the  optimization  of  am  actuation  system  for 


a  large  transport  aircraft,  while  General  Electric' e  work  (Air  Force  con¬ 
tract  AF  33(615) - 3587)  is  concerned  with  the  optimization  of  actuation  sys- 
tems  for  a  fighter  aircraft.  The  results  of  General  Electric's  work  have 
not  yet  been  published. 

The  Martin  Company  conducted  an  actuation  optimization  program  for  NASA 
entitled,  Optimization  of  Hydraulic  Thrust  Vector  Control  Systems  for  Launch 
Vehicles.  A  portion  of  this  work  is  reported  under  NASA  Accession  Num¬ 
ber  N66- 12611.  An  article  was  published  (Reference  3)  that  describes  the 
development  of  a  computer  program  for  related  system  cost  investigations. 
Costs  in  this  case  included  development  cost,  unit  cost,  weight  cost,  failures 
cost,  development  time  costs,  and  others.  Initially,  optimization  was 
obtained  by  observing  trends  in  computer  output  data,  rather  than  by  applica¬ 
tion  of  some  automatic  optimization  technique.  It  is  possible,  however,  that 
work  not  yet  reported  may  be  investigating  such  technique?!. 

AiResearcli  Manufacturing  Company,  Divisior  of  tne  Garrett  Corporation, 
conducted  an  cptiinizati'-.i  study  for  NASA  on  control  moment  gyros  (CMG's). 
The  CMG  design  was  optimized  with  respect  to  a  minimum  of  weight,  size, 
bearing  friction,  and  windage  losu.  Optimization  was  obtained  by  observing 
trends  in  computer  output  data.  A  report  of  this  work  is  given  in  Reference  4. 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


Based  upon  the  preceding  evidence,  the  following  conclusions  have  been 
drawn: 

1.  A  hydraulic  actuation  system  that  uses  current  design  techniques 
provides  a  near-optimum  design. 

2.  The  optimal  design  is  quite  dependent  upon  the  assumptions  and 
conditions  that  accompany  the  design  problem  so  the  individual  is 
still  very  much  a  part  of  the  design  process. 

3.  The  design  problem  used  in  this  study  was  beneficial  in  the  develop¬ 
ment  of  the  optimal  technique  but  the  ensuing  results  (because  of  the 
large  nominal  values  and  relative  insensitivity  of  the  cost  functions) 
did  net  constitute  an  adequate  base  upon  which  to  realize  the  full 
potential  of  the  design  technique. 

4.  It  is  possible  that  application  would  have  more  significance  to  high- 
performance  missiles  because  variations  in  actuation- system  per¬ 
formance  have  a  decided  effect  upon  fUght-control-system 
performance.  Also,  in  missile  systems,  it  would  be  possible  to 
conduct  actuation  power  system  tradeoffs. 

5.  It  appears  that  the  greatest  value  of  optimal  techniques  applied  to 
actuation- system  design  is  connected  with  comparison  of  different 
design  concepts  for  the  same  application.  The  study  of  various 
redundant  configurations  is  a  good  example. 

As  a  result,  it  is  recommended  that  additional  studies  be  conducted  to  obtain 
results  that  can  be  used  to  make  a  complete  evaluation.  This  work  should 
include: 

1,  Optimization  of  several  actuation  system  design  concepts  for  an 
actual  piloted  aircraft  flight  control  system. 

2.  Conducting  of  similar  studies  on  an  aerospace  vehicle  other  than  a 
piloted  aircraft. 

These  investigations  would  not  only  provide  optimal  designs  for  comparison 
with  designs  utilizing  current  methods,  but  would  provide  data  that  could  be 
used  to  study  the  problems  involved  in  the  direct  comparison  of  performance 
indexes.  Also,  the  work  involved  in  using  the  developed  technique  on  differ¬ 
ent  design  concepts  and  different  systems  can  be  determined. 
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APPENDIX  I 

STi.  ADY-STATE  MODEL  RELATIONSHIPS 


Each  of  the  sections  in  this  Appendix  contains  a  summary  of  th  parameters, 
output  equations,  and  calculations  used  to  generate  the  requireu  information 
for  the  following  (see  Section  II  for  further  reference): 


1. 

Pump, 

7. 

Valve  actuator  package. 

2. 

Accumulator. 

8. 

Fluid  temperature. 

3. 

Reservoir. 

9. 

System  elasticity. 

4. 

Tubing. 

10. 

Fluid  elasticity. 

5. 

Actuator. 

11. 

Actuator  cylinder  elasticity. 

6. 

Valve  assembly. 

12. 

Mechanical  linkage. 

1.  THE  PUMP 

Table  I  summarizes  the  parameter  relationships  associated  with  the  pump. 
Derivations  of  these  relationships  follow  the  summary  sheet. 

a.  Pump  Weight 

A  Vickers  slide  rule  can  be  used  to  determine  the  weight  for  various  pump 
sizes.  When  plotted,  the  pump  data  can  be  approximated  by  a  straight-line 
relationship  as  shown  by  Figure  27. 

Based  on  the  above  data,  pump  weight  as  a  function  of  pump  flow  can  be 
expressed  by 


PUW  =  7.5  +  0.615  (PJQ) 

The  system  has  four  tandem  valve -actuator  packages  and.  so,  has  eight 
valve-actuation  combinations.  Assuming  a  10%  internal  leakage  loss,  pump 
flow  af  a  function  of  valve  actuator  flow  can  be  expressed  by 

PUQ  •■=  i.l  (8 ) (Q) 


Combining  this  expression  with  the  previous  one,  results  in 


PUW 


7. 5  +  5.412  (Q) 


Table  1 
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D 

V  -  *0. +193. 7*0 

fciiASiury  j_ 

U 

%  -  4t5.+.0Hl*'**E8*Q 

WEIGHT  P 

Jf 

W  (PEES- 3000 y?0  -  7.343.412*0 
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Corporation,  Detroit,  ttlchlpan. 
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Figure  27.  Pump  Weight  Data 

Pump  weight  is  also  a  function  of  system  pressure.  Data  from  Vickers  indi¬ 
cate  an  approximate  weight  increase  of  5%  for  an  increase  in  pressure  from 
3,  000  to  4,  000  psi.  If  no  weight  reduction  is  assumed  for  pressures  less  than 
3,  000  psi  and  the  same  increase  is  assumed  for  every  1,  000  psi  above  3,  000 
psi,  then  the  pump  weight  expressed  as  a  function  of  valve  flow  and  pressure 
is 

PUW  =  [7.5  +  5.412  (Q)][l  +  5  x  10'5  (PRES  -  3,  000)1 


where 


(PRES  -  3,  000)  2  0 


(It  is  assumed  that  the  pressure  denoted  oy  PRES  is  a  nominal  pressure.  As 
a  result,  it  is  not  necessary  to  consider  the  pump  return  pressure  of  70  psi.) 

U  T"Y. -  r* 

v  .  JL  Utlip  VjUSI 

The  pump  co°t  of  the  various  sizes  and  for  an  arbitrary  quantity  can  be 
approximated  by  the  straight-line  relationship  shown  by  Figure  28. 
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Figure  28.  Pump  Cost  Data 

From  the  above  data,  pump  cost  as  a  ium.tion  of  valve  flow  can  be  expressed 
by 

PUC  *  1,400  +  11.3.2  (Q) 

c.  Pump  Size 

The  size  of  the  various  Vickers  pumps  expressed  as  cubic  inches  of  volume 
can  be  approximated  by  the  straight-line  relationship  shown  by  Figure  29. 

Based  on  the  data  shown  in  Figure  29,  pump  size  as  a  function  of  valve  flow 
can  be  expressed  by 

PUV  =  40  +  199.7  (O' 

d.  Pump  Reliability 

* 

The  pump  reliability  is  a  function  of  the  load  io  a  pump  hzi 

failure  rate.  If  horsepower  is  used  as  a  measure  of  the  load,  reliability  can 
be  correlated  with  flow  and  pressure.  Based  on  the  assumption  that  the  fail- 
ure  rate  is  reduced  about  10%  between  a  Vickers  3915  and  3911  pump  at  a 


pressure  of  3,000  psi,  pump  reliability  expressed  as  the  number  of  failures 
for  10&  operating  hours  can  be  approximated  by  the  straight-line  relationship 
shown  by  Figure  30. 


Based  on  the  data  shown  in  Figure  30,  pump  reliability  as  a  function  of  pres¬ 
sure  and  valve  flow  can  be  expressed  by 

PUR  =  465  +  0.01  11  (PRES)  (Q) 

2.  THE  ACCUMULATOR 

Table  II  summarizes  the  parameter  relationships  associated  with  the  accu¬ 
mulator.  Derivations  of  these  relationships  follow  the  summary  sheet. 

An  accumulator  is  used  to  damp  out  pressure  surges  caused  primarily  by 
sudden  opening  or  closing  of  the  metering  valve  and  to  smooth  out  the  ripples 
in  the  hydraulics  which  are  generated  by  the  pump. 

For  this  study,  accumulator  design  will  be  based  upon  the  additional  flow 

*■»•*•••  —  “  “  •  •*'»»!  fc.".  •«.  ••  MM  4  «•••«»  A.  —  ••  o  A J  ••  a.  Lot  -v  M  M  o  of 

o  wi  i ,y  i  cv^ui  i  cu  ww  c  v  wm  ojf  obcui  o  o  o  u«  c  ss  vm  u  - 

fied  value  after  a  sudden  valve  opening  and  for  the  period  of  time  the  pump 
takes  to  increase  its  output  flow  from  a  minimum  to  maximum  .alue.  Pres¬ 
sure  at  this  time  is  difficult  to  calculate  accurately  because  the  accumulator 


OUTPUT  PUMP  HORSEPOWER 


Figure  30.  Pump  Failure  Rate 

is  discharging,  which  reduces  pressure,  and  the  pump  flow  is  increasing, 
which  increases  the  pressure.  As  a  result,  the  minimum  design  pressure 
will  be  obtained  simply  by  averaging  initial  system  pressure  with  the  result¬ 
ing  accumulator  pressure  obtained  after  the  period  of  time  that  corresponds 
to  the  pump  time  constant.  Charging  or  compression  of  the  gas  in  the  accu¬ 
mulator  will  be  considered  an  isothermal  process  and  discharging  or  expan¬ 
sion  of  the  gas  will  be  an  adiabatic  process. 

The  pump  time  constant,  as  shown  by  Figure  31,  is  only  0.06  sec.  With  a 
control  surface  rate  requirement  of  40*/ see,  the  control  surface  will  travel 
less  than  1*  before  the  pump  will  be  able  to  put  out  100%  flew  when  starting 
from  5%  at  a  nominal  speed.  As  a  result,  the  minimum  acceptable  pressure 
when  using  flow  from  the  accumulator  will  be  assumed  to  be  2,400  psi.  This 
is  arbitrary  and  is  set  at  this  low  pressure  because  the  longest  it  can  last  is 
only  0.06  sec.  The  flow  required  from  the  accumulator  when  increasing  the 
demand  from  5%  to  100%  by  stepping  valve  position  is  shown  by  the  triangle 
on  the  previous  sketch  that  is  bounded  by  the  pump  flow  and  100%  flow  lines. 
Assuming  that  100%  flow  is  about  30  gpm,  the  accumulator  flow  is 
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Table  II 


ACCUMULATOR  SUMMARY  SHEET 


P#U£ 


/V/)A/£  :  .  accumulator _  5Yai6oL  _£ _ 


/ a/ Fu  ts :  JL _ L  J _ i _ OUPtlTS  *  c  c  w 

_  _A _ C _ C_ 

A  C  C  V 


..  '  ..tr.  L _ 

(A !£/£,#r  A 

C  C 

V 

-  3.-*(PRES-1000.)*.001 

Cost-  _a 

C  C 

c 

-  150. 

JI*£ 

C  C 

V 

-  90. 

FmtBujTy  ___* 

c  c 

Jt _ 

-  12.3 

// OFFS  '  Data  obtalnad  fro*  Bandix  Pacific  Dlvlalon,  Bandix  Aviation  Corporation, 

North  Hollywood,  California. 


Figure  31.  Pump  Response  Time 


ACCQ  =  Area  of  triangle  ABC 

=  3.  47  in."* 

With  a  5-in.  spherical  accumulator  precharged  to  2,000  psi  and  an 
volume  of  60  in.^,  an  isothermal  filling  to  3,  000  psi  will  result  in  a 
volume  of 


P1V1 


P2  v2 


AAA  (  /  Al 


\OU|  =  J,VWJ  \Vp) 


initial  gas 
.  new  gas 


If  3.4?  in.  of  fluid  15  used,  the  resulting  pressure  following  an  adiabatic 
expansion  will  be 


P  V 
2  2 


1.4 


P,  = 


P,  = 


P  V 
3  3 


1.4 


3’000  (iot) 

2, 670  psi 


1.4 


The  minimum  calculated  pressure  would  be,  therefore, 


P 

min. 


3,  000  +  2,  670 

- 2 - 


=  2, 835  psi 

which  is  well  above  the  minimum  allowable  pressure.  This  particular  accu- 
r.vilator  also  can  be  used  with  systems  of  larger  flow  capacity.  Working 
backwards,  starting  with  a  minimum  allowable  pressure  of  2,400  psi  and 
maintaining  a  2,000  psi  precharge,  100%  flow  can  be  as  high  as  153  gpm. 

This  type  of  accumulator  is  quite  flexible  and  is,  therefore,  independent  of 
system  flow.  As  a  result,  its  cost,  size,  and  reliability  will  be  assumed 
constant  and  will  have  the  following  values: 


Cost  -  $ 

ACCC  =  150. 

Size  -  in . ^ 

ACCV  =  90. 

Reliability  -  failures 

ACCR  =  12.  3 

per  10&  operating  hours 

The  weight,  however,  is  a  function  of  pressure  and  can  be  approximated  by 
the  straight-line  relationship  shown  by  Figure  32. 

Based  on  the  above  data,  accumulator  weight  as  a  function  of  supply  pressure 
can  be  expressed  by 


ACCW  =  3+  (PRES  -  1,  000)(10'3). 
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Table  111  summarizes  the  parameter  relationships  associated  with  the  reser¬ 
voir.  Derivations  of  these  relationships  follow  the  summary  sheet. 

a.  Reservoir  Weight 

The  reservoir  size  is  dependent  on  the  quantity  of  fluid  in  the  system. 
Because  the  tubing  contains  most  of  the  fluid,  reservoir  design  is  primarily 
a  function  of  tube  size.  Weight,  therefore,  can  be  calculated  by  a  standard 
Douglas  calculation  used  for  preliminary  sizing,  as  follows: 

1 .  Assume  a  tubing  fluid  weight  (TUFLW)  and  then  add  a  number  of 
factors . 


?.  Add  thermal  expansion  fsc*or  =  5.5%  (TTTFT.Wi. 

3.  Add  leakage  factor  =5%  (TUFLW). 

4.  Add  accumulator  fluid  weight  =  ACCFLW. 


TABLE  HI 

RESERVOIR  SUMMARY  SHEET 


PAt£ 
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SIZE 

^ - 

_JL_ 

R 

E 

W 

-  1.8*.215*TUFUI 

C 
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* 
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R 
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R 

E 

L 
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R 

E 

D 
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/Sorssi 


Dct*  obtained  fro*  Aircraft  Division,  Douglas  Aircraft  Co., 
uong  Bench,  Callfc.-nla. 


Based  on  the  above  data,  reservoir  weight  as  a  function  of  tubing  fluid  weight 
can  be  expressed  by 

REW  =  1 , 8  +  0.219  (TUFLW) 

Because  return  pressure  will  not  vary  appreciably  and  because  the  major 
portion  of  the  reservoir  is  concerned  with  return  pressure,  it  will  be 
assumed  that  the  reservoir  constraints  are  not  functions  of  the  supply 
pressure. 

b.  Reservoir  Cost 

The  variation  in  cost  will  depend  on  the  change  in  the  amount  of  material 
required  because  the  same  configuration  will  be  used  in  each  case.  Even 
though  this  change  is  minor,  the  relationship  shown  in  Figure  34  will  be  used 
to  determine  retervoir  cost. 

Based  on  the  data  shown  in  Figure  34,  reservoir  cost  as  a  function  of  tubing 
fluid  weight  can  be  expressed  by: 

REC  =  245  +  0.91  (TUFLW) 


5  200 


TUBING  FLUID  WEIGHT  (LB) 


Figure  34.  Reservoir  Cost  Data 


c.  Reservoir  Sise 


Be  reuse  of  the  odd  shape  of  the  reservoir  (a  large  diameter  at  only  one  end 
and  a  quite  small  diameter  at  other  end),  the  large  diameter  and  overall 
length  will  be  considered  instead  of  the  total  volume.  Assuming  that  a  diam¬ 
eter  and  stroke  combined  change  occurs  when  its  capacity  changes,  the 
reduction  in  these  dimensions  for  a  change  in  capacity  will  be: 

Letting 

V  =  KD2S 

where 

V  =  reservoir  fluid  capacity. 

D  *  large  reservoir  diameter. 

S  =  reservoir  piston  stroke. 

K  =  proportionately  constant. 

Then  a  50%  reduction  in  capacity  will  result  in 


Using  standard  reservoir  design  procedures  and  a  straight-line  relationship, 
the  calculated  variations  in  length  and  diameter  can  be  shown  by  Figure  35. 

From  the  above  data,  reservoir  length  and  diameter  as  a  function  of  tubing 
fluid  weight  can  be  expressed  by 

REL  =  9  +  0.18  (TUFLW) 

RED  =  1  +  0.04  (TUFLW) 
d.  Reservoir  Reliability 

Reservoir  reliability  will  be  assumed  constant  because  the  same  configura¬ 
tion  will  be  used  in  each  case.  DC-8  reliability  data  indicate  a  failure  rate 
for  boots  trap -type  reservoirs  to  be  on  the  order  of  16  failures  per  10”  oper¬ 
ating  hours.  Therefore: 


RER  =  16 
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0 


Figure  35.  Reservoir  Size  Data 


4 ,  TUBING 

Table  IV  summarizes  the  parameter  relationships  associated  with  the  tubing. 
Derivations  of  these  relationships  follow  the  summary  sheet. 

a.  Tubing  Weight 

Standard  Douglas  tube  sizes  used  for  a  3,  000-psi  system  are  shown  in 
Table  V. 

The  corresponding  plumbing  weight  in  lb/ft  which  includes  fitting,  clamps, 
and  fluid  can  be  approximated  by  the  straight-line  relationship  shown  in 
Figure  36. 
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lane  laad,  California. 


•  » 


Aircraft  Co, 
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Table  V 

STANDARD  PLUMEING  SIZES 


Wall  thickness 

Tube  size 

Steel 

Aluminum 

(in.) 

(supply  system) 

(return  system) 

0.035 
0.042 
.049 
.065 


_  1-0 


STEEL 


M  Qi 

TUBWG  (OO-M.) 


Fiffire  38.  Plunking  Weight  Data 


Based  on  the  above  data,  plumbing  weight  in  lb/ft  of  length  for  a  3,000-psi 
system  as  a  function  of  the  tubing's  outside  diameter  can  be  expressed  as 
follows: 

For  steel 

TUWSF  =  0,311  +  1.6  (DS  -  0.5) 

For  aluminum 

TUWALF  =  0.162  +  0.695  (DAL,  -  0.5) 

where 

DS  and  DAL  2  0.5 

As  the  plumbing  is  divided  into  two  categories,  the  total  plumbing  weight  will 
be: 

Total  tubing  weight  =  No.  1  steel  +  No.  2  steel  +  No.  1  aluminum 

+  No.  2  aluminum 

TUW  =  TUWS  1  +  TUWS  2  +  TUWAL  1  +  TITW  AL  2 

where 

TUWS  1  =  (0.311  +  1.6  (DS  i  -  0.5)1  TULS  1. 

TUWS  2  =  [0.311  +  0-16  (DS  2  -  0.5)1  TULS  2. 

TUWAL  1  =  [0.162  +  0.695  (DAL  1  -  0.5)]  TULAL  1. 

TUWAL  2  =  [0.162  +  0.695  (DAL  2  -  0.5)1  TULAL  2. 

To  maintain  about  the  same  stress  levels  in  the  tubing  at  various  system 
pressures,  wall  thicknesses  will  change,  and  this,  in  turn,  will  affect  the 
weight.  Currently,  all  the  required  combinations  of  standard  wall  thick¬ 
nesses  and  tube  diameters  for  the  higher  pressure  systems  are  not  available. 
However,  it  will  be  assumed  that  they  would  become  available  if  the  need 
existed  in  quantities  to  make  their  manufacture  worthwhile.  The  diameter 
and  wall  thickness  combinations  which  could  bo  used  for  the  steel  pressure 
lines  are  shown  in  Table  VI. 


77 


Table  VI 

PLUMBING  SIZES  AS  A  FUNCTION  OF  SYSTEM  PRESSURE 


Tube  size 
(in.) 

3,  000 

Pressures  (psi) 

4,  000  5,  000 

Wall  thickness 

6,  000 

1/2 

0.035 

0.049 

0.058 

0.065 

5/8 

0.042 

0.038 

0.065 

0.083 

3/4 

0.049 

0.065 

0.083 

0.095 

1 

0.065 

0.083 

0.095 

0.120 

A  weight  calculation  indicates  an  increase  of  about  12%  at  all  tube  sises  for 
each  1,000-psi  increase  in  pressure.  As  a  result,  the  plumbing  weight 
expressed  as  a  function  of  tubing  diameter  and  system  pressure  is 

TUW  =  JTUWS  1  +  TUWS  21  (1  +  0.00012  {PRES  -  3,  000)1 
+  TUWAL  1  +  TUWAL  2 


where 


(PRES  -  3, 000)  £  0  . 

Because  the  return  pressures  will  not  appreciably  change  regardless  of  the 
supply  pressure,  it  will  be  assumed  that  the  return  system  plumbing  weight 
will  not  change  as  a  function  of  pressure. 

Because  the  reservoir  design  is  a  function  of  the  fluid  weight  in  the  tubing, 
the  weight  of  the  fluid  will  be  calculated  separately. 

Using  the  following  expression  between  the  tubing  inside  and  outside  diameter 
for  a  3,  000-psi  system. 

TUID  ■  0.430  +  0.88  (D  -  0.5) 


where 


D  2  0. 5  , 


7t 


i 


the  fluid  weight  per  foot  of  length  expressed  as  a  function  of  the  tubing  outside 
diameter  is 


TUFLWF  =  0.293  {0.430  +  0.88  (D  -  0.5)12 

Fluid  specific  gravity  was  assumed  constant  at  a  value  of  0.86.  It  was  also 
assumed  that  the  effect  of  the  difference  of  wall  thickness  between  the  steel 
and  aluminum  tubing  for  a  3,  000-psi  system  is  negligible,  so  that  the  above 
expression  would  apply  to  both  types  of  tubing. 

Separating  each  tube  category  results  in  the  following  expressions: 

TUFLW  =  TFWS  1  +  TFWS  2  +  TFWAL  1  +  TFWAL  2 


where 

TFWS  1 
TFWS  2 
TFWAL  1 
TFWAL  2 


0.293.(TULS  1)  {0.430  +  0.88  (DS  1  -  0.5))2 
0.293  (TULS  2)  {0.430  +  0.88  (DS  2  -  0.5))2 
0.293  (TULAL  1)  {0.430  +  0.88  (DAL  1  -  0.5))2 
0.293  (TULAL  2)  {0.430  +  0.88  (DAL  2  -  0.5)]2 


Use  of  the  higher-pressure  systems  with  the  thicker-walled  tubing  results  in 
less  fluid  weight.  Referring  to  Table  VI,  the  reduction  in  inside  diameter  for 
each  1,  000-psi  increase  in  pressure  regardless  of  tube  size,  averages  out  to 
be  4.9%.  As  a  result,  (TUFLW)  with  a  T'-essure  term  added  will  be 


TUFLW  =  (TFWS  1  +  TFWS  2]  (1  -  2.4  x  10'9  (PRES  -  3,  000)2] 
+  TFWAL  1  +  TFWAL  2 


where 


(PRES  -  3,  000)  2  0  . 


b.  Tubing  Cost 

Costs  for  tubing,  fittings,  clamps,  and  an  arbitrary  quantity  can  be  approxi¬ 
mated  by  the  straight-line  relationship  as  shown  in  Figure  37. 

Based  on  the  foregoing  data,  plumbing  cost  in  dollars  per  foot  of  length  for  a 
3,  000-psi  system  as  a  function  of  the  tubing's  outside  diameter  can  be 
expressed  as  follows: 

For  steel 


TUCSF  =  1.3  (DS) 


TU6M6  (00  -  H.) 


Fifare  37.  Plwbiaf  Cost  Data 

For  ziuminum 


TUCALF  =  0.222  (DAL.) 

As  «  result  of  the  division  of  (dumbing  lengths,  the  total  cost  will  be 


TUC  =  TUCS  1  4-  TUCS  2  +  TUCAL  1  +  TUCAL  2 


where 

TUCS  1  =  1.3  (DS  1)  (TULS  1) 


TUCS  2  =  1 .3  (DS  2)  (TULS  2) 

TUCAL  1  =  0.222  (DAL  l)  (TULAL  1) 
TUCAL  2  =  0.222  (DAL  2)  (TULAL  2) 


i 

} 


Calculations  indicate  about  a  13%  increase  In  cost  for  each  1,000-pai  increase 
in  pressure  above  3,  000  psi.  As  a  result,  the  plumbing  cost  expressed  as  a 
function  of  tubing  diameter  and  system  pressure  is 

TUC  =  [TUCS  1  +  TUCS  2]  (l  +  0.00013  (PRES  -  3,  000)} 

+  TUCAL  1  +  TUCAL  2 


where 


(PRES  -  3, 000)  2  0  . 


c .  Tubing  Size 

Tubing  volume  (and  the  space  it  occupies)  normally  does  not  affect  the  corre¬ 
sponding  vehicle  envelope  requirements  for  this  type  of  aircraft;  therefore, 
these  factors  will  not  be  considered  in  this  study. 

d.  Tubing  Reliability 

Tubing  reliability  will  be  considered  constant  regardless  of  tube  size  because 
the  thicker  wall  tubing  is  used  with  the  higher  pressures.  Thus,  for  each 
case,  about  the  same  maximum  stress  level  is  retained.  A  failure  rate 
obtained  from  a  FARADA  Handbook  for  a  horizontal  stabilizer  hydraulic  tub¬ 
ing  assembly  is  44  failures  per  1G#*  operating  hours.  Therefore, 

TUR  =  44 


e.  Tubing  Pressures  Lass 

The  tubing  pressure  losses  and  fluid  temperatures  are  interrelated.  There¬ 
fore,  to  calculate  these  losses,  a  fluid  temperature  will  be  assumed.  Because 
this  is  a  simplified  hydraulic  system  and  to  reduce  the  complexity  of  the 
pressure  loss  calculations,  the  fluid  temperature  will  be  assumed  to  be  a 
constant.  Based  upon  the  hydraulic  system  design  capabilities  generated  dur¬ 
ing  the  C-5A  proposal  effort,  a  fluid  temperature  of  120°F  will  be  used. 

With  120°F  fluid  and  flows  in  the  low-turbulent  region,  the  pressure  loss  can 
be  approximated  by  a  straight-line  relationship  on  a  log-log  plot,  as  shown 
by  Figure  38. 

Based  on  these  data,  tubing  pressure  loss  per  foot  of  length  as  a  function  of 
tubing  flow  and  tubing  outside  diameter  for  a  3,  000-pai  system  can  be 
expressed  by 


TUDPF  =  (TUQ)1'75  [0.00155  (D)"4,93] 

It  was  assumed  that  the  effect  of  the  difference  in  Wall  thickness  between  the 
ste^l  and  aluminum  tubing  for  a  3,  000-psi  system  is  negligible,  so  that  the 
above  expression  would  a*  ply  to  both  types  of  tubing. 


PRESSURE -FLO*  RELATIONSHIP  (AP/  1.7$) 


Adding  10%  of  (TUDPF)  for  pressure  losses  in  plumbing  bends,  fittings,  and 
other  small  components  not  considered  in  this  study,  and  separating  each 
tube  category,  results  in  the  following  expressions  for  line  loss  as  a  function 
of  valve  flow,  line  length,  and  line  diameter. 

TUDP  =  DPS  1  +  DPS  2  +  DPAL  1  +  DPAL  2 

where 

DPS  1  =  0.001705  iTULS  1)  (8.8Q)1,75  (DS  1)“4'93 

DPS  2  =  0.001705  (TULS  2)  (4.4Q)1’75  (DS  2)"4,93 

DPAL  1  =  0.001705  (TULAL  1)  (8.8Q)1'75  (DAL  if4,93 

DPAL  2  =  0.001705  (TULAL  2)  (4. 4Q)1-75  (DAL  2)“4'93 

Use  of  the  higher-pressure  systems  with  the  thicker -walled  tubing  results  in 
smaller  flow  areas  and  more  pressure  loss.  Referring  to  Table  III,  the 
reduction  in  inside  diameter  for  each  1,000-psi  increase  in  pressure  regard¬ 
less  of  tube  size  averages  out  to  be  4.9%.  The  corresponding  increase  in 
pressure  loss  is  about  30%.  As  a  result,  (TUDP)  with  a  pressure  term  added 
will  be: 


TUDP  =  [DPS  1  +  DPS  2]  [1  +  0.0003  (PRES  -  3,000)] 
+  DPAL  1  +  DPAL  2 


where 


(PRES  -  3, 000)  >  0  . 
f.  Tubing  Size  Limitation 

Because  a  relatively  high  fluid  temperature  was  used  for  the  tubing  pressure 
loss  calculations,  selection  ov  a  low  limit  on  tube  size  was  necessary  to 
ensure  that  the  system  would  have  sufficient  actuator  AP  to  perform  satisfac¬ 
torily  at  lower  temperatures.  For  this  application,  a  low  temperature  of 
-40°F  has  been  selected,  and  satisfactory  performance  at  -40°F  has  been 
defined  as  being  able  to  obtain  about  25%  of  maximum  control  surface  rate 
with  zero  hinge  moment.  This  means  that  the  tubing  pressure  loss  at  -40°F 
i  annot  exceed  the  difference  between  the  supply  pressure  and  the  required 

vy  V P>  AP  Tn  an  artnal  annlirafirtw  frK  o  trfnati rvr»  eirofom  should 

— *■  ■  -  — - — rx - - - *  - - -  /  -  * - -  r - * - - - — 

be  checked  at  various  fluid  temperatures  after  the  tubing  sizes  have  been 
selected  to  ensure  that  the  system  meets  the  different  satisfactory  perfor¬ 
mance  requirements  whatever  they  maybe. 

Using  a  laminar  flow  viscosity  correction  factor  for  a  nominal  pressure  of 
3,  000  psi,  pressure  loss  data  ran  be  approximated  by  the  straight-line 


relationship  plotted  on  log-log  paper  and  shown  by  Figure  39.  Baaed  on  theae 
data,  the  pressure  loss  at  «40°F/ft  of  length  for  a  3.  OOQ-psl  system  as  a 
function  of  tubing  flow  and  outside  diameter  can  be  expressed  by: 

TUDPF  =  (TUQ)  (0.306)  (D)"3, 98 

It  was  assumed  that  the  effect  of  the  difference  of  wall  thickness  between  the 
steel  and  aluminum  tubing  for  a  3,  000-psi  system  is  negligible  so  that  the 
above  expression  would  apply  to  both  types  of  tubing. 

Adding  10%  of  (TUDPF)  for  pressure  losses  in  plumbing  bends,  fittings,  and 
other  small  components  not  considered  in  this  study  and  separating  each  tube 
category  results  in  the  following  expressions  for  line  loss  as  a  function  of 
valve  flow,  line  length,  and  line  diameter. 

TUDP4  =  DP314  +  BPS24  +  DPAL14  +  DPAL24 


where 


DPS14 

=  2.69  (TULS1)  (Q)  (DSl)-3’98 

DPS24 

=  1.35  (TULS2)  (Q)  (DS2)-3'98 

DPAL14 

=  2.69  (TULAL1)  (Q)  (DALI)”3,98 

DPAL24 

=  1 .35  (TULAL2)  (Q)  (DAL2)"3,98 

Use  of  the  higher-pressure  systems  with  the  thicker -walled  tubing  results  i» 
more  pressure  losses  because  of  smaller  flow  areas  and  a  larger  viscosity 
correction  factor.  Referring  to  Table  VI,  the  reduction  in  the  inside  diame¬ 
ter  for  each  1, 000-psi  increase  in  pressure  regardless  of  tube  size,  aver¬ 
ages  out  to  be  4.9%.  The  corresponding  increase  to  pressure  loss  is  about 
20%  because  of  the  tubing  and  also  20%  because  of  the  viscosity  correction. 
As  a  result,  (TUDP4)  with  a  pressure  term  added  will  be 

TUDP4  r  [DPS14  +  DPS24J  [1  +  0.0002  (PRES  -  3. 000)1* 

+  D  PALI  4  +  DPAL24 


where 


(PRES  -  3,000)  2:  0 


5.  ACTUATOR 


Table  VII  summarise!  the  parameter  relationships  associated  with  the 
actuator,  f  erivations  of  these  relationships  follow  the  summary  sheet. 

a.  Actuator  Weight 

A  schematic  of  the  tandem  actuator  unit  is  thc?m  by  Figure  40 .  Sample 
weight  calculations  bated  upon  this  schematic  and  plotted  as  a  function  of 
torque  arm  length  with  the  supply  pressure  as  a  parameter  indicate  a  rela¬ 
tionship  as  shown  by  Figure  41 .  Application  of  a  hyperbolic  relationship  to 
these  curves  with  a  minimum  weight  correction  factor  results  In  the  follow- 
in  expression: 


ACTW  =  HMT  (10"5)  10.0625  i  +  4  {RNOM 


+  140  (PNOM)"0' 242 


To  obtain  the  above  expression,  the  ratio  of  the  actuator  AP  (ACTDP)  to  the 
supply  pressure  (PRES)  was  approximated  by  a  value  of  two-thirds.  This 
ratio  was  used  because  it  approximates  the  pressure  where  maximum  power 
transfer  occurs.  Because  the  system  flow-pressure  relationships  are  inters 
related,  an  actuator  AP  approximation  was  necessary  to  calculate  the  actua¬ 
tor  area. 

In  the  above  expression,  actuator  design  is  based  upon  a  constant  torqua 
design.  During  the  optimization  process,  however,  the  actuator  weight  lor 
variable  torque  design  will  have  to  be  determined  because  arbitrarily 
selected  valuer  of  supply  pressure,  actuator  area,  and  torque  arm  length 
will  be  used.  Corresponding  torque  correction  factors  must  be  added  to  this 
expression.  Witn  *»<*o  of  the  same  type  of  calculations,  the  actuator  area  fac¬ 
tor  can  be/4pproximated  by 

+3.3  (AREA-ANOM)  (10-5) 

The  supply  pressure  factor  can  be  approximated  by 

+0.95  (PRES-PNOM)  (10*8) 

The  torque  arm  factor  can  be  approximated  by 


+  1.2  (R-RNOM)  (10-5) 
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ACTUATOR  SUMMARY  SHEET 
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Figure  40.  Actuator  Schematic 

Incorporation  of  the  (HM)  notation 

HM  =  (1/2)  HMT 
results  in  the  final  expression 


ACTW 


1  +  4  ^RNOM  -  5  6 (PNOM] 
+  140  (PNOM)'0,242  +  3.3  (AREA-ANOM) 

+  0.95  (PRES-PNOM)  (10"3)  +  1.2  (R-RNOM)] 


=  2  HM(10‘5) 
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b.  Actuator  Cost 

The  actuator  it  an  integral  part  of  the  valve  actuator  package  so  its  coat  will 
be  included  in  the  total  assembly. 

c.  Actuator  Size 

Because  the  valve  assembly  more  or  less  surrounds  the  actuator,  the  actua- 
tor's  area  is  not  an  important  dimension.  The  length  is  predominant,  how* 
ever,  and  will  be  considered.  The  actuator  length  is  primarily  a  function  of 
the  stroke.  Other  factors,  such  as  method  of  attachment,  cylinder  ends, 
piston  thickness,  etc.,  influence  actuator  length;  however,  for  a  particular 
configuration  they  will  be  fairly  constant.  As  a  result,  actuator  length  will 
be  defined  in  terms  of  a  constant  plus  a  stroke  function. 

Assuming  a  simple  small  angle  relationship,  actuator  stroke  as  a  function  of 
control  surface  deflection  can  be  expressed  by 

ACTS  =  DELMX  (R)  /57. 3 


Because  the  actuator  is  a  tandem  unit,  the  piston  rod  must  be  at  least  three 
times  as  long  at  the  calculated  stroke  (two  cylinders  +  one  stroke).  For  the 
Boeing  727  elevator-control  package,  it  has  a  length  factor  of  about  13  in. 
over  and  above  that  contributed  by  the  stroke.  As  a  result,  the  overall  actu¬ 
ator  length  can  be  represented  by: 

ACTL  =  0.0524  (DELMX)  (R)  +  13 
d.  Actuator  Reliability 

As  with  the  other  components,  published  reliability  data  are  very  general, 
and  ao  it  is  impossible  to  use  it  to  obtain  a  correlation  between  failure  rates 
and  actuator  design  parameters.  To  determine  whether  actual  failure  reports 
would  provids  more  useful  information,  it  was  decided  to  survey  the  ATA 
(Air  Transport  Association)  failure  summary  reporta  and  determine  whether  a 
correlation  axlste,  for  example,  betw  t  actuator  leakage  failures  and  actu¬ 
ator  stroke  because  stroke  i  a  parameter  being  used  in  this  study.  ATA 
reports  covering  the  period  of  l  January  1961  through  31  December  1964  for 
the  DC-8  totaling  1,271,201  hours  of  aircraft  operation  were  reviewed.  The 
results  are  tabulated  in  Table  VIII. 
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Table  VIII 

DC-8  ACTUATOR  LEAKAGE  FAILURE  DATA 


Function 

No.  per  plane 

Nominal  stroke  a 

No.  of  failurea 

Main  landing  gear  uplatch 

2 

0.6 

3 

Aileron 

2 

4.1 

0 

Tab  lockout 

2 

4.2 

1 

Flap  outboard 

2 

4.3 

0 

Flap  midwing 

2 

5.2 

0 

Flap  inboard 

2 

7.1 

2 

Main  landing  gear  trim 

2 

7.4 

1 

Noae  wheel  ateering 

2 

7.9 

1 

Rudder 

1 

7.9 

0 

It  ia  not  evident  from  the  above  data  that  there  ia  a  correlation  between  leak¬ 
age  failurea  and  atroke  if  one  doea  exiat.  A  factor  that  complicatea  the  aitua- 
tion  ia  the  incorporation  of  deaign  changea  during  the  4-year  period  becauae 
it  waa  impoaaible  to  determine  the  change  effectively  from  the  available  data. 
It  waa  noted  during  the  aurvey,  however,  that  when  a  failure  waa  not  diag- 
noaed  ra  an  iaolated  incident  and  a  deaign  change  occurred,  thia  change 
uaually  corrected  the  aituation  and  no  more  failurea  of  thia  kind  happened 
again. 

Aa  a  reault,  becauae  the  actuator  ia  an  integral  part  of  the  valve  actuator 
package,  ita  reliability  will  be  included  in  the  reliability  of  the  total  aaaembly. 

6.  VALVE 

Table  DC  aummarizea  the  parameter  relationahipa  aaaociated  with  the  valve. 
Derivationa  of  theae  relationahipa  follow  the  aummary  aheet. 

a.  Valve  Aaaembly  Weight 

The  valve  aaaembly  containa  a  number  of  componenta.  Theae  include  the 
following: 

1.  Tranafer  valve. 

2.  Autopilot  actuator . 


TABLE  IX 

VALVE  SUMMARY  SHEET 
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FACTATt _ _ 


5yAfSoL  -1_  —  —  — 


/fi/fUTS:  - 

V  _A.JL.JJ _ 


QU.FU.TS 


v _ a _ f _ JL 

V  A  y  .  c  - 

_* _ A _ l - £- 
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8—  ACT _ 

-  16. _ 


Data  obtained  tr*»  B.rtae  Corporation,  Irvin.,  California 
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3.  Metering  valve. 

4.  Tab  lockout  valve. 

5.  Filters. 

6.  Solenoids. 

7.  Bypass  valves. 

8.  Input,  feedback,  and  summing  leakages. 

9.  Various  manifolds . 

Because  this  unit  is  similar  to  the  Boeing  727  elevator -valve  assembly,  and, 
because  it  contains  so  many  different  functions,  its  weight  will  be  expressed 
as  a  function  of  valve  flow  and  based  upon  the  weight  of  the  Boeing  727  unit. 
The  Boeing  unit  weighs  about  35  lb  and  has  a  maximum  valve  flow  of  about 
1  gpm,  so,  if  10%  weight  is  added  for  each  gpm  above  one,  the  valve  assem¬ 
bly  weight  for  a  3,  000-psi  system  can  be  expressed  by 

VASW3  =  35  +  (Q  -  1)  (0.10)  (35) 


For  systems  with  pressures  other  than  3,  000  psi,  a  nominal  change  in  weight 
of  5%  for  each  1,  000-psi  change  in  system  pressure  will  be  used.  Adding  this 
pressure  function  to  the  previous  expression  results  in 

VASW  *  (31.5  +  3. 5Q)  [1  +  0.00005  (PRES-3,  000)J 

b.  Valve  Assembly  Cost  and  Reliability 

Because  the  valve  assembly  is  an  integral  part  of  the  valve  actuator  package, 
its  cost  and  reliability  will  be  included  in  the  total  unit. 

c.  Valve  Assembly  Size 

The  valve  assembly  size  is  quite  flexible  because  its  combination  of  valves, 
manifolds,  and  so  forth,  can  be  patched  together  in  many  different  possible 
wrya.  As  a  result,  it  can  be  made  to  fit  most  envelopes  and  so  the  volume  or 
size  of  the  unit  will  not  be  considered. 

7.  VALVE  ACTUATOR 

'fable  X  summarizes  the  parameter  relationships  associated  with  the  valve 
actuator.  Derivations  of  these  relationships  follow  the  summary  sheet. 
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TABLE  X 

VALVE  ACTUATOR  SUMMARY  SHEET 


/V/)A/£ :  -  valvs _  5y/“£oL  J2 _ i _ I _ 

AIIML! 


/A/fars:  _L  _L  _L  _JL _ Oumrs  jL_i.-L.jL 

Q  _  _  _ 


nurpur  jrm/Ar/OA/s _ 

£//.#}'  XX-L  Jl _ •  OM  *  «-«—  " 

Cost"  _ i~*» _ 

J/*£  JLJ.JLJ.___  _ _ 

foot eury _ Jm.w  _ 


tfO/ y?S\  IMiaklilMtl  tnm  ScrtM  Ctr^riltw,  trviM,  Ciltfinti 


a.  Valve  Actuator  Package  Weight 

'.his  ur.it  is  a  combination  of  the  valve  assembly  and  actuator  and  so  its 
weight  can  be  expressed  by 

VAPW  =  ACTW  +  VA3W 

b.  Valve  Actuator  Package  Cost  and  Reliability 

Cost  and  reliability  will  remain  about  the  same  for  the  various  systems  that 
are  studied  because  the  same  configuration  will  be  used  in  each  case.  Data 
obtained  from  Bertea  Products  using  the  Boeing  727  unit  as  a  guide  indicate 
that  for  an  arbitrary  quantity,  the  cost  in  dollars  is 

VAPC  =  7,  500 

The  un  .  reliability  based  upon  data  obtained  from  various  commercial  air¬ 
lines  for  the  Boeing  727  unit  over  a  period  of  20  months  and  a  total  service 
time  of  188,441  flight  hours  indicates  the  number  of  failures  per  10°  hours 
to  be 


/ 

VAPR  =  16 


c .  Valve  Actuator  Package  Size 

The  predominant  size  factor  in  this  assembly  is  the  actuator  length.  As  men¬ 
tioned  for  the  valve  assembly,  the  overall  size  is  not  considered  because  the 
design  is  very  flexible  and  can  usually  be  made  to  fit  most  envelopes. 

d.  Valve  Actuator  Package  Pressure  Distribution 

System  pressure  can  be  divided  into  the  following  three  parts:  (1)  plumbing 
AP,  (2)  actuator  AP,  and  (3)  valve  AP,  as  shown  by  Figure  31.  The  flow- 
actuator  rate  relationship  is  a  function  of  the  actuator  area  which,  in  turn,  is 
a  function  of  the  actuator  AP.  As  shown  by  Figure  42,  however,  the  available 
actuator  AP  depends  upon  the  flow.  As  a  result,  an  assumption  must  be  made, 
values  calculated,  the  assumption  compared  to  the  calculated  value,  and  the 
process  repeated  if  the  comparison  reveals  too  large  an  error  A  method  that 
can  be  used  to  obtain  a  good  first  assumption  is  a 3  followb: 

At  maximum  rate 


PRES  =  TUDPM  *  VASDPM 


DESIGN  REQUIREMENT  =  0  HM  AT  MAXIMUM  RATE 


Pressure  Distribution 


because  there  isn't  any  actuator  AP.  If  it  is  assumed  that  the  plumbing  AP  is 
a  function  of  instead  of  then  at  one-half  rate  the  actuator  AP  can 

be  defined  as 


ACTDP1  =  PRES  -  (1/4)  TUDPM  -  (1/4)  VASDPM 
=  PRES  -  1/4  [TUDPM  +  VASDPM] 

=  3/4  PRES 


With  a  value  for  actuator  AP,  actuator  area  at  the  one-half  rate  design  condi¬ 
tion  can  be  expressed  by 


AREA  =  HMD/  ACTDP1  (R) 


The  corresponding  valve  flow  in  gallons  per  minute  can  be  expressed  by 

Q  =  0.00453  (DET.D)  (AREA)  (R) 


From  this  information,  a  new  value  for  the  actuator  AP  (ACTDP)  at  the  one- 
half  *  rate  design  condition  can  be  calculated 

ACTDP  =  PRES  -  0.25  (PRES  -  TUDPM)  -  TUDP2 

To  determine  whether  an  iteration  is  required,  the  following  ground  rule  was 
established: 

if  (ACTDP1)- (ACTDP)  <  100,  then  use  ACTDP  as  calculated, 
if  (ACTDP1)- (ACTDP)  >  100,  then  substitute  (ACTDP1  +  ACTDP)/2 
for  the  old  value  of  ACTDP1  and  recalculate  ACTDP. 

8.  FLUID  TEMPERATURE 

A  heat  exchanger  was  not  included  in  this  simplified  actuation  system  because 
a  cursory  fluid  temperature  check  indicated  that  the  fluid  temperature  would 
not  exceed  the  nominal  specified  limit  of  160*^.  The  following  calculations 
are  a  sample  of  the  type  used  to  check  the  fluid  temperature.  When  used  on  a 
DC-8,  temperatures  obtained  utilizing  these  calculations  compared  favorably 
with  actual  test  data.  The  following  assumptions  will  be  used: 

1.  Overall  pump  efficiency  and  all  system  pressure  looses  with  the 
exception  of  actuator  AP  go  into  heating  the  fluid. 

2,  Nominal  flow  of  one-third  maximum  pump  flow. 


3.  Actuator  AP  of  three-fourths  system  pressure. 


4. 


All  of  the  heat  lost  will  be  by  convection  with  a  surface  coefficient 
(C)  of  2  Btu/hr  ft2  °F. 

5.  Surface  area  will  consist  of  tubing  area  plus  10%  for  other 
components. 

6.  A  standard  sea  level  ambient  day  temperature  IT)  of  68°F  will  ere 
used  because  ground  operation  will  give  the  higher  fluid 
temperatures . 

With  a  hinge  moment  of  36,  000  in.  -lb,  a  control  surface  rate  of  40°/ sec,  and 
a  system  pressure  of  3,  000  psi,  the  nominal  system  flow  (QI'T)  will  be 

nM  (0.00453)(DELD)(HMD)(l/3)(8.8) 

=  (3/4)(PRES) 

(0 . 00453) (40)(36,  000)(l/3)(8. 8) 

=  (3/ 4} (3,  000) 


=  8.5  gpm 


The  hydraulic  system  output  power  (HPO)  based  upon  the  actuator  flP  will  be 

_  13/4)(PRES)1QN) 

"  1714 

(3/ 4)(3,  000)(8. 5) 

1714 

=  11.15  hp 


The  input  power  (HPI)  with  an  overall  pump  efficiency  of  92%  will  be 


HPI  = 


1714  (0.1 

■ 

=  16. 15  hp 


Therefore,  the  heat  generated  (HIN)  will  be 


Rt;i  /  V>  r 

HIN  =  (HPI-HPO)  ^44 1  ~  fig— 

=  (16.15  -  il .  15)(254?) 

=  12,730  Btu/hr. 


The  tubing  surface  area  (TSA)  will  be  calculated  on  the  basis  of  using  a  5/8- 
in.  -diam  for  the  No.  1  tubing  and  1/2-in.  -diam  for  the  No.  2  tubing.  Aa  a 
result,  the  total  surface  area  (SA)  will  be 

SA  =  1.1  (TSA) 

=  ((0.625)  (370)  +  (0.5)  (160)} 

=  89.5  ft2 

l 

The  fluid  temperature  obtained  with  12,  73u  Btu/hr  heat  input,  89.5  ft2  of  sur¬ 
face  area,  and  68° F  ambient  will  be 

H1N  =  (C)  (SA)  (FLT  -  T) 

12,730  =  2(89.5)  (FLT  -  68) 

FLT  =  1 39°  F 

9.  SYSTEM  ELASTICITY 

Table  XI  summarizes  the  parameter  relationships  associated  with  the  sys¬ 
tem  elasticity.  Derivations  of  these  relationships  follow  the  summary  sheet. 

The  actuation  system  resonant  frequency  is  a  function  of  the  system  elastic¬ 
ity,  so  relationships  between  it  and  the  system  design  parameters  must  be 
obtained.  In  this  study,  those  items  which  will  be  used  to  describe  system 
elasticity  are  fluid  elasticity  and  the  elasticity  of  the  actuator  cylinders. 

They  are  like  springs  in  series  and  the  effective  elasticity  of  the  system  can 
be  expressed  by 


1  1  .  1 
SE  s  flE  +  EyE 


10.  FLUID  ELASTICITY 

Table  XII  summarizes  the  parameter  relationships  associated  with  the  fluid 
elasticity.  Deri  vations  of  these  relationships  follow  the  summary  sheet. 

The  fluid-bulk  modulus  (or  elasticity)  is  a  function  of  the  fluid  temperature 
and  pressure.  Data  obtained  by  R.  L.  Feeler  and  J,  Green  (Reference  5)  on 
the  adiabatic  bulk  mnHnln;  of  MIL-H-56C6A  fluid,  covering  a  temperature 
range  of  50  s  to  200  #F,  can  be  approximated  by  Figure  43.  Based  on  these 
data,  the  fluid  bulk  modulus  as  a  function  of  the  fluid  temperature  and  system 
pressure  can  be  expressed  by 

FLE  =  -610  (FLT)  +  138500  f(2/3)(PRES)]0, 1082 
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TEMPERATURE-SULK  MODULUS  RELATIONSHIP!  FLE  +  .61  (FLT)  (10*)  1  10"*) 


In  this  expression,  it  was  assumed  that  the  actuator  cylinder  pressure  at 
steady  state  and  for  small  cyclic  inputs  can  be  approximated  by  a  value  of 
two-thirds  of  the  system  pressure. 

11.  ACTUATOR  CYLINDER  ELASTICITY 

Table  XIII  summarizes  the  parameter  relationships  associated  with  the 
actuator  cylinder  elasticity.  Derivations  of  these  relationships  follow  the 
summary  sheet. 

Under  static  conditions,  expansion  of  the  actuator  cylinder  tends  to  lower 
cylinder  pressure  so  it  will  be  considered  in  determining  the  system  elastic¬ 
ity.  Cylinder  expansion  is  a  function  of  the  following: 

1.  Cylinder  pressure. 

2.  Piston  area. 

3.  Cylinder  thickness . 

4.  Piston  rod  diameter. 

5.  Material  modulus  of  elasticity  (E). 

6.  Poissons  ratio. 

Assuming  that  the  last  three  factors  are  constant  and  using  stress-strain 
relationships  for  a  thick-walled  steel  cylinder,  cylinder  elasticity  can  be 
approximated  by  the  straight-line  relationship  shown  by  Figure  44.  Based  on 
these  data,  the  actuator  elasticity  for  a  steel  cylinder  as  a  function  of  system 
pressure  and  actuator  area  can  be  expressed  by 


=  0.0019  (AREA)-0 • 41  (PRES)"0,88 
12.  MECHANICAL  LINKAGE 

Table  XIV summarizes  the  parameter  relationships  associated  with  the 
mechanical  linkage.  Derivations  of  the  relationships  follow  the  summary 
sheet. 

The  mechanical  linkage  includes  all  of  those  parts  used  in  the  transmission  of 
a  mechanical  algncl  front  the  cockpit  to  the  valve  actuator  package  for  each 
elevator  control  surface.  This  includes  cockpit  control  columns,  cable  tension 
regulators,  torque  overside  mitt  l.o.i  feel  mechanisms,  snd  miscellaneous 
cabling  and  linkages.  This  mechanical  control  link  will  be  considered  as  a 
complete  assembly  with  weight  being  the  only  item  of  interest.  Because  none 
cf  the  design  parameters  being  considered  influences  the  design  of  this  con¬ 
trol  link,  its  weight  will  be  considered  constant  &t  a  v*Ju-  cl 

MLW  =  108  lb 


m 


TABLE  XIII 

ACTUATOR  CYLINDER  ELASTICITY  SUMMARY  SHEET 
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TABLE  XIV 

MECHANICAL  LINKAGE  SUMMARY  SHEET 


APPENDIX  II 


DYNAMIC  MODEL  RELATIONSHIPS1 


A  hydraulic  actuation  system  operates  in  a  nonlinear  way  and  so,  logically, 
it  should  be  represented  hy  a  nonlinear  dynamic  model.  Such  a  model  is 
shown  in  Figure  45.  The  model  includes  the  essential  characteristics  of  an 
actual  hydraulic  flow-control  system.  Those  refinements,  however,  such  as 
valve  lap,  valve  leakage,  and  actuator  torque-arm  elasticity,  which  identify 
a  particular  system,  have  not  been  included  for  simplicity.  The  absence  of 
these  refinements  influences  model  performance  primarily  in  the  resonant 
area;  however,  as  shown  in  Figure  46,  a  portion  of  the  characteristic  reso¬ 
nant  peak  is  still  present.  Data  on  a  more  complete  model  can  be  found  in 
Reference  6. 

When  designing  a  control  system,  some  form  of  electronic  compensation 
normally  is  used  to  counteract  the  adverse  effect  of  the  resonant  peak  on 
system  performance  and  to  more  or  less  linearize  system  response.  In  this 
design,  a  linearized  response  and  a  linearized  model  were  desired  to  facili¬ 
tate  use  of  the  ISE  criterion  but,  as  mentioned  previously,  the  electronic 
portion  of  the  control  system  was  not  included  and  so  could  not  be  used  for 
this  purpose.  As  a  result,  to  retain  the  essential  hydraulic  system  charac¬ 
teristics,  and  still  have  a  fairly  linear  response,  it  was  decided  to  add 
dynamic  pressure  feedback  to  the  actuation  system.  For  simplicity,  the 
effect  of  the  dynamic  pressure  feedback  will  be  shown  only  on  the  final  linear¬ 
ized  model.  Because  the  dynamic  pressure  feedback  will  not  affect  the  equa- 
f5ons  derived  in  this  section,  its  effect  can  be  added  to  the  final  form  by 
application  of  the  superposition  principle. 

Referring  to  the  nonlinear  dynamic  model  (Figure  45),  displacement  of  the 
valve  s pool,  X,  is  the  input  to  the  system,  and  the  actuator  piston  motion, 

Y  ,  is  the  output.  From  the  block  diagram,  flow  into  cylinder  No.  1  is 

Qcl,9)  ■  KlKzJPs  -  pcl«  X<s> 


If  Xn  is  defined  a*  the  nominal  spool  position,  Pcn  the  nominal  cylinder  pres¬ 
sure,  and  small  perturbations  about  these  values  are  considered,  flow  is 
given  as 


3Q  .(s)  9Q_,(s) 

*Qcl<S)  =  3P^*Pcl<3>+ 


2K1K2 


K  -  ^i^)-17  Vpci<R) + kik2(p.  -  pcnr-<.) 


1/2 


1  The  symbols  used  in  Appendix  II  have  been  simplified  as  an  aid  to  the 
development  of  the  required  functions.  All  symbols  used,  however,  have 
been  defined  in  the  appendix. 
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WHERE 

X  -  VALVE  SPOOL  MOTION 

K1  -VALVE  AREA -STROKE  RELATIONSHIP 

K2  -  FLOW  CONSTANT 

Av  -  VALVE  METERING  AREA 

Qg  -  VALVE  FLOW 

Vc  -  UNCOMPRESSED  FLUID  VOLUME  SUBJECT  TO  CYLINDER 

Ps,  P„-  SUPPLY  AND  RETURN  PRESSURE 

V,  -  ACTUATOR  VOLUME  SUBJECT  TO  CYLINDER  PRESSURE 

Pe  -  CYLINDER  PRESSURE 

/3  -  SYSTEM  ELASTICITY 

Ap  -  ACTUATOR  PiSTON  AREA 

I  -  LOAD  INERTIA 

Yp  -  ACTUATOR  OUTPUT 

I  -LOAD  OUTPUT 

Ypu  -  MAXIMUM  ACTUATOR  OUTPUT 

A  P  -  ACTUATOR  PISTON  DIFFERENTIAL  PRESSURE 

R  -  ACTUATOR  MOMENT  OR  T  ORQUE  ARM 

!  -  OPWTOfi  d/dt 


Figure  45.  Nonlinear  Dynamic  Model 


AMPLITUDE 


FREQUENCY  (CPS) 


The  nomine!  conditions  are 


Combining  this  with  thr  representation  of  the  oynamic  pressure  feedback 
given  by  Moog  (Reference  7)  the  block  diagram  shown  in  Figure  47  describee 
the  linearised  system. 

To  verify  accuracy  of  the  linearised  modei,  a  frequency  response  of  the 
iinaar  model  »»•  compared  with  the  frequency  response  of  the  nonlinear 
model  taken  at  differ  ont  input  levels.  The  results  for  a  10%  input  is  shown 
in  Figure  46.  As  is  readily  seen  frotn  these  data,  the  frequency  response 
curves  of  ^a  nonlinear  model  with  dynamic 'pressure  feedback  and  the  linear 
model  compare  well.  From  additional  data  taken,  but  not  shown,  good 
comparison  was  found  for  input  levels  up  to  30%. 


X'.  -  COMMAND  INPUT 

X  -  VALVE  SPOOL  MOTION 

K.  -  AMPLIFIER  AND  TORQUE  MOTOR  GAIN 

Kq  -  FLOW  CONSTANT 

Kp  -  PRESSURE  FEEDBACK  GAM 

r  -  TNE  CONSTANT  OF  PRESSURE  FEEDBACK 

Kj  -  FEEDBACK  GAIN 

\  -  EQUIVALENT  HYCRAUUC  SPRING 

A  P  -  ACTUATOR  PISTON  DIFFERENTIAL  PRESSURE 

Ap  -  ACTUATOR  PISTON  AREA 

I  -  LOAD  INERTIA 

Yp  -  ACTUATOR  OUTPUT 

8  -  LOAD  VELOCITY 

R  -  ACTUATOR  MOMENT  OR  TORQUE  ARfc 

S  -  OPERATOR  d/dt 


I 

i 


AMPLITUDE 


FREQUEACY  (CPS) 


APPENDIX  III 


NO.  1  COMPUTER  PROGRAM 


1.  DESCRIPTION  OF  THE  FIXED-GRID  COMPUTER  PROGRAM 

The  computer  program  was  written  with  a  series  of  subroutines  for  the 
evaluation  of  the  cost  functions.  This  was  done  for  ease  of  understanding  the 
completed  program  and  for  convenience  in  making  cc-rections  and  additions. 
The  entire  program  is  written  in  Fortran  IV  language.  The  main  program 
basically  sets  up  the  grid  over  which  the  performance  index  will  be  evaluated 
and  calculates  the  performance  index.  There  are  six  subroutines:  one  each 
for  the  five  cost  functions  and  one  for  the  evaluation  of  the  nominal  system. 
The  fixed  grid  is  used  to  evaluate  in  the  following  manner: 

a.  The  program  fixes  a  pressure,  moment  arm,  and  the  tube  sizes, 
and  then  evaluates  the  minimum  area  that  will  meet  the  hinge- 
moment  requirements.  The  performance  index  and  cost  functions 
will  be  evaluated  for  these  parameter  values. 

b.  Then  holding  all  of  the  parameters  fixed,  the  area  will  be  increased 
by  0.  5  sq  in.  and  the  performance  criterion  re-evaluated.  This  will 
be  repeated  three  times. 

c.  Then  the  pressure  will  be  increased  by  adding  a  fixed,  predetermined 
amount  (250  psi  for  this  example),  then  the  minimum  area  will  again 
be  calculated,  and  the  performance  index  evaluated.  The  area  will 
then  be  incremented  as  indicated  in  step  b. 

d.  After  the  entire  range  of  pressure  has  been  investigated  at  the 
initial  moment  arm  and  tube  diameters,  the  moment  arm  will  be 
decreased  by  a  fixed  constant  (0.  5  in.  for  this  example).  Then  the 
entire  range  of  pressures  and  the  corresponding  areas  will  be 
re-evaluated.  This  will  continue  until  all  of  the  moment-arm  values 
to  be  investigated  have  been  used, 

e.  After  this  has  been  accomplished  for  the  initial  values  of  the  tubes, 
the  entire  process  must  be  repeated  again  for  each  possible  combi¬ 
nation  of  the  tube  sizes.  This  will  cover  the  parameter  space. 

The  logic  for  each  subroutine  is  similar.  The  flow  diagram  for  the  weight 
subroutine  shown  by  Figure  49  is  typical  of  this  logic. 


TOTAL  >  MW/ 
.AND  <  MAX. . 


*  RETURN  TO 
kMAtN  PROCRAI 


Figure  49.  weight  Subroutine  Logic 


2.  COMPUTER  VARIABLES  NOT  PREVIOUSLY  DEFINED 


Computer  variables  not  previously  defined  and  u»ed  in  tfce  No.  1  listing 
include  the  following: 


TOTW 

total  weight. 

TNOMW 

total  nominal  weight. 

TOTC 

total  cost. 

TNOMC 

total  nominal  cost. 

TOTDYP 

total  dynamic  performance. 

TOTV 

total  volume. 

TNOMV 

total  nominal  volume. 

TOTR 

total  reliability. 

TNOMR 

total  nominal  reliability. 

VEQ 

equivalent  volume. 

VEQN 

nominal  equivalent  volume. 

REDN 

nominal  reservoir  diameter. 

RELN 

nominal  reser^^r  length. 

ACTLN 

nominal  actuator  length. 

WR 

desired  bandwidth  of  actuation  system. 

LISTING  FOR  NO.  I  COMPUTER  PROGRAM 


-IT  FOR  MAIN 


OPTIMAL  DESIGN  OF  AN  ACTUATING  SYSTEM 
DIMENSION  PARAMt 31 )«  FJ5t6 ) 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 


/OAT At  At  38  I 

/NOM J /  tnomw  ,tnomc*tnomv*tnomr*tnomoy  *t nomen 

/O UP/  TOTW»TOTC*TOTV*ToTR.TOTDVP«TOTENVtACTW 
/FLG/  IEPElG*IERFLG*!DIA.NNEND.NN 


mainoooo 

MAIN0003 

MAIN0006 

MAIN0009 

MAIN0012 

MAIN001S 


/  !NTM/  DS1  *  OS  2*  DAL  1  *DAL2  *  Q»  PRES»AREA»R»ACTOP*TUF>.W«ANORM  MAIN0018 


EQUIVALENCE  t  At  26 > *TUu«L i ) . t  A 1 2? > «TULAL2 ) ♦ ! At  28 ) * TULS1 ) ♦ 
1  <  AI2«>)  »TULS2) 

WRITE  16*633) 

READ  t  3«300)  tAtI)«I*l*38) 

WRITE  1 6*330)  tAtI)*I«l*38) 

READ  15*501)  NFND*  MEND 
WRITE  16*551)  N'NO.  MEND 
WRITE  1 6*555) 

DO  10  II*-  2*31 

10  PARAMt 1 1 1  «  0.0 
IEPELG  •  0 
IERELG  «  0 
NNEND  -  5 
IENO  *  1 
JFMR  ■  1 
KEND  *  1 
LENO  •  1 

IE t  A(7).NF.  0.0)  LENO  •  2 
IE t  At  9) «NE.  0.0)  KENO  •  2 
IF t  Aflll.NE.  0.0)  JFNO  •  2 
IF  t  A 1 13 ) *NE.  0.0)  IENO  •  2 
PARAMt 1)  r  AI3I 
0013  11-2*9 

13  PARAMt 1 1 )  ■  PARAMt I I-l )  ♦  A13) 

PARAMt 1)*5. 5 
PARAMt 10)  *  A 1 2 ) 

DOlB  11*2*17 

18  PARAMt 1 1+9)  ■  PARAMt 1 1 ♦S 3  ♦  At*) 

PARAMt 10) »  A  1 1 ) 

00  200  I*  1  * IENO 


00  200  I* 
DO  200  J* 
DO  200  Km 
DO  2C0  L* 
DO  200  M* 
DO  200  N* 
DO  199  NN* 


J*  1 tJEND 
Km  1 *KEND 
L*  l.LENO 
M*  l.MEND 
N*  ltNEND 
NN*  1.  NNEND 
PARAMt  N*9) 


MAIN0019 
MAINC020 
MAIN0023 
MAIN002* 
MAIN002T 
MAIN0030 
MA IN0033 
MA 1X003* 
MAIN0036 
MASN0QJ9 
►MIN0062 
MAINOOAS 
MAIN0069 
MAIN0072 
MA1N0073 
MAIN0078 
MAIK0081 
MAIN008* 
MAIN0087 
MAIN0090 
MAIN0093 
MA IN0096 

MAIN0102 


R  *  PARAM(M) 

AREA  ■  AREA  ♦  .3 
DAI  1  *  AtL  +  5) 

DAL2  •  AtJ49) 

DS1  •  AiK+7) 

DS2  •  AtJ+11) 

CAuL  NOMSYS 

IE t  ID IA  *EQ.  0)  GO  TO  23 
213  •  PRES  -  3000, 

!EtZ23.LT.O.O)  ZZ3  *0.0 

ZZl  r  3366* t  TULSl*tDSl)**t-3.98)*tl.+.0002*ZZ3)**2 
U)**t-  J.38 ) ) 

ZZ2  «  •  3366* t  TULS2*tDS2)**t-3,9l»*n.*.0002*ZZ3)**2 
12)**t-3.98)) 

23  TUDP*  •  0*  t  2.2*221  ♦  ZZ2) 

IE  t  PRES  .GF.  TUPP*  )  GO  TO  Id 
1U 


MAIN011* 
MAIN011T 
MAINC120 
MA1K0123 
MAIN0126 
MAIN0129 
MAIN0112 
MA I NO 139 
MAIN0162 
MAIN0163 
MAIN0168 
MA  I  NO 1 7 1 
MAIN017* 
MAIN017? 

U  A  flirt  9  *  ft 
nw  m-  •  m •* 

MAIN01S6 


TULALl*fOAL 
TULAL2*  t  OAL 


MAIN013* 
MA  II  9136 


WRITE  (6*3941  I  »J»K*L*M»N»NN*TUDP4 
IF  ( IPPFLG  .tO.  0 )  GO  TO  210 

go  to  m 

30  CAUL  WEIGHT 

IFf  IERPLG  «EQ«  Ql  SO  TO  33 

32  PIDX  «  0.0 

IERFLG  •  0 

IF  < IFPFLG  ,E0,  0)  GO  TO  210 
GO  TO  115 

33  CALL  COST 

IF  {  IERFLG  .ME.  01  GO  TO  32 
CALL  DYNPER 

IF  1  IERFLG  .ME.  01  GO  TO  32 
CALL  VOLUME 

IF  f  ( ERF! G  «NE.  0)  GO  TO  32 
CALL  RELIA8 

IF  (  IERFLG  .ME.  0  I  GO  TO  32 
IF  UFPFLG.NE.  0)  GO  TO  110 
PARAMC 1  !■> ( 3  I 
PARAM( 10) «  A  (2  1 
IFPFLG  ■  1 
110  PIDX  •  0.0 

FJSdl  ■  1.0  -  TOTW/(2«0*TNOMVn 
FJS<2»  «  1.0  -  TOTC/(2.0*TNOMC1 
FJSdl  ■  1.  -  TOTOYP 


MAIMOSSft 
HA  INI 136 
MAIN0197 
MAIN01J* 
NAINOUi 
MAIN© 144 
MA I NO 14? 
MAIN014G 
MAIN0190 
MAIN0130 
MAIN0132 
MAIMU3? 
MAIN2132 


MA1N0139 


K41N0140 
MAIN0193 
MA  I  NO 134 
MA INI 134 


FJS(4»-  1 .O-TOTV/ ( 2.0*TNOMV  J 
FJ5<51  •  1.0  -  TOTR/  (2.0  *  TNOMR 1 
00  112  II  »  1.3 

IF  HFJSmt.LT.  0.01.0R.1FJS1III.GT.1.01 1  GO  TO  115 

112  CONTINUE 

00  1U  II  •  1.3 

114  P10X  •  PIOX  ♦  A!II+31>*FJS(III 

113  WRITE  (6*5321  I .J.K.L.M. M.NN.PIOX.TOTW. TOTC. TOTOYP. TOTV.ToTR. 
1ACTW 

199  CONTINUE 
NNEND  *  * 

200  continue 

210  STOP 

300  FORMAT  (  7F10.3  1 
501  FORMAT  (  2(3X.  12 1  1 

330  FORMAT  (  1H0.  ATX.'JHINPUr  DATA////13H  P  NOM  -  F10.3.11H  P 


MAIN&194 

MAIN7134 

MA1N8194 

MAIN9134 

MAIN0133 

MAIN013* 


MAIN01I0 
MAIN0K9 
MAIN0192 
MAIN0193 
MAIN0196 
MAIN0201 
MINNA IN0204 

Is  •  F10.?.10H  R  NOM  ■  F10.3.10H  P  INC  •  F10.3.HM  R  INC  »  MAIN0207 
2F10.3//13N  0AL11  »  F10.3.11H  DAL  12  -  F10.3.10H  0S11  •  MAIN0210 

3F1O.3.10H  OS  12  ■  F10.3.11H  0AL21  ■  F10.3.11H  0AL22  »  F10. 3//MA INQ21S 

413H  DS21  »  F  10.3.  llH  DS22  «•  F10.3.10H  W  MAX  *  F10„3.10H  WMAINQ216 

3  MIN  «  F10. 3. I 1H  C  MAX  ■  F10.3.11H  C  MIN  •  F10.3//J3H  PUV  MAXMAIN0219 


6  -  F10. 3*1 1H  RED  MAX  >  F10.3.10H  R  MAX  ■  F10.3.10H  R  MIN  -  MAIN0222 
TF10.3.11H  CYP  MAX  -  F10.3.11H  OYP  MIN  ■  F10.3//13N  ENV  MAX  •  MAIN0223 
8F10.3.11H  ENV  MIN  «  F10.3.10M  TULAL1  •  F10.3.I0H  TULAL2  •  F10.3.  MAIN022* 
91 lH  TULS1  -  F10.3.11H  TULS2  "  F10.3//13H  DELTA  MAX  ■  F10.3.  MAIN0231 
111H  MM  «  F10.3.10H  WTF1  ■  F10.3.10H  C0STF2  •  F10.3.11H  OYPFMA IN0232 


23  ■  F10.3//13H  V0LF4  ■  F10.3.11H  RELF5  •  F10.4.10H  ENVF6  MA1N12S2 

3“  F10.3.11M  REL  MAX  ■  F10t?i  MAIN1232 

331  FORMAT  (  JHO*  11H  P  STATES  ■  IA.3X.10HR  STATES  •  14 1  MAIN023A 

352  FORMAT  (1H  .  6H  CASE  511. 12. 1 1 • 3X.7HP. ! •  ■  F7.2.3X. 

19MWT  ■  FS.1.3X.7HC0ST  «  F7.0.3X.6HDVP  ■  F6.2.SX.6HV0L  •  F6.2.3X.  MAIN0240 
26HREL  «  F7.1.3X.7HACTW  FT. 21  MAIN1240 

333  FORMAT  (1H1.  AOX.  2AH0PTIMAL  ACTUATING  SYSTEM  1  MAIN0241 

93A  FORMAT  (  1M  .  6H  CASE  1 1 » II . 1 1 . 1 1. 1 1 . 1 2 . 1 1 .3X.  6HTUBP4  -  KAIN0242 

1  F1C.31  KAIN02A9 

393  FORMAT! 1H1.34X.36HCASE  »  DS2.DAL2.0$liDALi.R. PRES. ARE*  1  MAIN124S 
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“IT  *0*  HOHSYS 

SUBftOU?  ■  iic 


ZZJ  -  on?;,?-'1  60  TO  3U’0,LS<‘-.PSEi,"RF.  -  .  XOMSOO 

28! :  «r  ,AC70P>T— "0*H  «bs 

OALlO  •  0,0  NOMSOOji 

0AL20  •  o.o  HOHSOQl; 

GO  To  s  NOMSOOIJ 

3  fF  'NN  .eo.  ..  NOMSOOil 

Ir ft NNEND.ro  l,\  00  TO  5  WOMS0012 

GO  To  J5  *  »ANO,  (Nn.£q  j  r.  "OMS0012 

5  }JM  .  HM/(PRCS#  Dt  M  G0  70  3  ^50013 

ANORw  «  Apr i  R  ^  NOMSOQ12 

Ir  '0AL1  .Fn\..  NOMSOOIS 

0&L10  a  04i  j  GO  TO  7  *QM$0018 

GO  To  10  #0M$0021 

•'  t  OSl.FQ  ncn  W0MS0121 

a,?o*s,«« ’ ” Go  ro  •  ssss 

OAL2ODiL0Ai°tDAL2O>  60  TO  9  *25oo|s 

GO  To  !2  12  N0MS002 6 

9  IF  (  osp  m  rr  ^®HS00>7 

0S20  »  0S7^* °3?0 )  GO  TO  m  NOMS0028 

1n  GO  To  I?  NOMSOoll 

10  I0IA  a  j  *OMS0030 

126  ■  pope  NO«SC* 

'J'zz*.tT.0;0?' ™:  :«e« 

*21  ■  ,00i*ni*5Z6  *  0.0  *&*SO 

,<lo"o",,,S”  ZULS1*(Q5i  j##  "0MS!> 

ZZj>  «  nn  laj>1*(DAL 

*  tuuet.,0  SS»- 

'OOP  .  i;?r,?'**l.T5  ,u,-*>-Z*{0*L 

tuopj  _  T. 

*GTDP  .  ppZj#  ZZ3  ♦  ill  #f*  ZZ|*ZZ,  WOMS004 

lF  <irPFLG  Np  *£5#<  P«ES  -  N0MS004 

■  <  .75 ;2fes°»  00  TO  25  °P  *  *  7W>P*  "OMSOoT. 

?r  'ZZ4.LMOO  «  *  ACTDP)  '•OMSOOs 

PRES  .  pres1?0;®1  GO  TO  20  N0MS005< 

,  GO  TO  5  5  +  ZZ*'  2.0  WOMS0037 

20  WRITE  (a,  NOMS003* 

23  RETURN  °'  PR£S>AREa,  a  rt°M5006/ 

90  FORMAT (  iho  .u  °*  TUDP,TtJ0P2,ACrnp  S°*S0063 

1  »h  tu^ho:  «h i  pres  ,  Fl0#,  ,ACr°p  "WisooiJ 

PN0  ’l0H  TUDP2  .  *  RIO.,.  ,u  „ 

*"”•  ’H  ACT0P'-  no,?/0’3*  SoSSSi! 

N0MS0076 
NOMSOOfl] 


•'vnjugi 

w&asooi 

nomsooj 


XOtisoo' 

klAlia  —  _ 


-IT  FOR  WEIGHT 

SUBROUTINE  WEIGHT 

COMMON  /DATA/  A(13l .WMX.WMN*  DUH( iO) .TULALl *TULAL2*TULS1*TUIS2» 

1  DELMX.  HM 

COMMON  /FLG/  IFPFLG* IERFLG* 10 1 A 

COMMON  / INTM/  OSi *052* DALI »DAL2 *Q* PRES, AREA  »R *AC TOP .TUFLW.ANORM 
COMMON  /OUP/  TOTW.OUMA(5) .MCTW 
COMMON  /NOMJ/  TNOMW 
IF  ( iFPFLG.NE.Ol  GO  TO  5 

112  «=ARS( SORT (0,06?5*(1 .+4  .*(  5 . 0-56.*  300l<  .**( -0, 280 1 1  **?  )  5 1*140, • 
13000. **(-0.242) 

ACCW  <*  3*  (PRFS  -  1000,1*0,001 
5  ZZ3  -  PRES  -  3000.0 
ZZ4  •  ZZ3 

IF f  ZZ4.LT.  CrO)  ZZ4  »  0.0 

PUW  »  (7.5  ♦  5.412*01*  ( 1.0+. 00005*  ZZ41 

VASW  ■  (31.5  ♦  3.5*0* » <  1,0  ♦  ,00005*ZZ3) 

IF  (  IOIA.EO.O)  GO  TO  10 

ZZ5  ■  (1.6*  0S1  -  ~»89)*TULS1  ♦  (1.6*  OS 2  -  .489)*TULS2 

ZZ6  ■  (  .695*0AL1  -  . 1855 ) *TULAL 1  ♦  (  ,695*DAL2  -  . 1855) *TULAL2 

TFWSl  *  , 293*1 ULS 1  * ( ,4  30  *  .*8* (DS1-. 51 1**2 

TFWS2-.293*TULS2*(.430+.38*(DS2-.5) 1**2 

TFwAL2«.293*TULAL2*( .430+. 88*( 0AL2-. 511**2 


WEIGOOOO 

WEIGOOO* 

WEIGOOOC 

WEIG9006 

WEIG0009 

WEIG0012 

WEIG0015 

WESG0C16 


WEIG0021 

WEIG0022 

WEIG0023 

WEIG0024 

WEIG0027 

WEIG0030 

WEIG0033 

WEIG0016 


TFWAL1»o293*TULAL1*(  .430+. 8 8* (PALI-.  5  1 1**2 

TUFLW  ■  (TFWSl  +TFWS2 1  * ( 1 •  ♦  2.4*10.** i -9 , 1 *( Z24 1 »*2 1  +  TFWAL1  ♦  T 
1FWAL2 

REW  »  1.8  ♦  .219*TUFUW 

10  ACTW».96* (ZZ2+3. 3* ( AREA-3.20 >+.95*{PRES“3000, 1/1 000. +1. 20* i R-5.0 ) ) 
TUW  -  ZZ5*  (1.0  ♦  ,00012*  ZZ41  ♦  lib 

XMLW  -  108. 

TOTW-PUW  +ACCW+  4.0*VASW  +  TUW  +  REW  ♦  4.0*ACTW  ♦  XMLW 

IF  ( 1FPFLG.E0.  0)  TNOMW  «  TOTW 

IF  i  (  TOTW. LE.WMX1. AND.  (  TOTW, GE.WMNP  GO  TO  15 

IFRFLG  ■  1 

WRITE  (6*25) 

15  RETURN 

25  FORMAT  (1H  .  lOX.  20HWEIGHT  OUT  OF  LIMITS  > 

END 


WEI60051 

WEIG0055 


WEIG0060 

WEIG0066 

WEIG0069 

WEIG0072 

WEIG0075 

WEIG0078 

WEIG0081 


-It  FOR  COST 

SUBROUTINE  COST 

COMMON  /DATA/  AU5)  .CMX.CMN.DUM16 )  .TULMJ  .TULAL2 .TULSl »TUiS2 
COMMON  /FL3/  IFPFLG. 1E3FLG. 1DIA 

COMMON  /INTM/  DSl  .D52.DAU  .DAL2 .0.PRES. AREA.* .ACTDP »TU?l.a 
COMMON  /OUP/  TCTW.TOTC 
COMMON  /NOMI/  TNOMW.TMOMC 
IF  UFPFLG.NF*  0)  60  TO  5 
ACCC  •  150.0 
VAPC  ■  7500,0 
5  PUC  «  1*00,0  +  123.2*  0 
IF  (  IDIA  »EQ.  0)  CO  TO  10 
IOI A  «  0 

REC  »  2*3.0  ♦  .91*  :FLW 

221  «  1.3  *«  0S1*TULS1  ♦  DS2»TULS2) 

222  •  .222  *<  0AL1*TULAL1  +  DAL2*TULAL2 I 

1C  223  «  PRES  -  3000.0 

IF  (  223  ,LT.  OeO>  Z23  ■  0,0 

TUC  -  2Z1»!  1.0  ♦  .00013*2231  ♦  222 

TOTC  «  ACCC  ♦  VAPC  ♦  PUC  <*■  REC  ♦  TUC 

IF  (  IFPFL6  .EQ.  Oi  TNOMC  •  TOTC 

IF  I  I  TOTC.LE.CMX > .AND. < TOTC. GE ,CMN> }  60  TO  13 

IERFLG  ■  1 

WRITE  16.25) 

15  RETURN 

25  FORMAT  <1H  ,10X,  1AHC0ST  OUT  OF  LIMITS  ) 

END 


COST 0000 
COSY 0003 
COST 0306 
COST 0009 
COSTOOU 
C0ST0019 
COSTOOiS 
C0ST0021 
COST 002* 
COSTCO?? 
COSTOCJO 
COST 0033 
C0ST0036 
COSTOOS9 
COSTOOA2 
C0S?00*5 
COST 00*9 
COST 0051 
COST 005* 
COST0057 
C0STC060 
COST0063 
COST 0366 
COST 0069 
COST0072 
C0ST0075 


-FT  FOR  OYNPER 

SUBROUTINE  OYNPfR 
COMMON  /DATA/  DUM <211 »DYPMX * DYPMN 
COMMON  /F LG/  IFPFLG.  IERFLG 
COMMON  / 1 NTM/  DUM2  <3)»PRES*ARtA*R 
COMMON  /OOP/  DUM3UI. TOTDYP 
COMMON  /KOMI/  DUMA  ( t> )  .TNOMDY 
<*  10. 

CY£  ■  AREA  **0.41*PRES**0.88/0()0019 
FIT  ■  120. 

FLE  ■  -610. ''FLT+138500„»(2.*PF.ES/3*)**.  1062 
SE  “  FLE  »  CYE/IFLF  +  CYE ) 

SIGM  ■  A0./57.3 
AYE  -  420. 

WN  ■  SORT  (2.*  BE  *ARFA*R / 1  SI GM  *  AYE)) 

TOTDYP  -  l.-EKPI“.25MWR-0.1*WN) ) 

IF (TOT OYP.LT .0.0)  TOTDYP  ■  0.0 
IFUOTDYP.GT.l.Ol  TOTDYP  «  1.0 
12  IF  ( IFPFLG. EO.  0)  TNOMDY  «  TOTDYP 

IF  (  (TOTDYP  .LE.  DYPMX ) .AND. ( TOTDYP  .GF,  DYPMN))  GO  TO  15 
IERFLG  »  1 
WRITE  (4.20) 

15  RETURN 

20  FORMAT  (IF1  .  10X.  23HDYN.  PER.  OUT  Or  LIMITS) 

END 


OYMPOOOO 

DYNP0003 

DYNP0006 

DYNP0012 

DYNP0015 


DYNP0024 

DYNP0027 

DYNP0030 

DYNP0033 

DYNP0036 

DYNP003R 

DYNP0042 
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-n  FOR  VOLUME 

SUBROUTINE  VOLUME  VOwUOOQO 

COMMON  /DATA/  A«17»*  PUVMXf  REDMX*  DUMCIOI.  DELMX*  DUH1  (7)  *RELMX  VO.UOOOJ 


COMMON  /FLG/  IFPFLG*  IERFL6 

V0UJ0004 

COMMON  /INTM/  DUM2<M,  Q.PRES»AREA.R*ACTDP»TUFLW 

VOwUOOOf 

COMMON  /OUP/  TOTW*TOTC*TOTV 

VOwUOOlZ 

COMMON  /NOM!/  TNOMW*TNOMC .TNOMV 

VOwUOOl# 

IF  (  IFPFLG  *NE.  0)  00  TO  9 

VOwUOOlS 

Z2 1  •  #092396  •  DELMX 

VOwUOOZl 

9 

PUV  •  40*0  ♦  199.7  *  0 

V0wU002A 

IF  (  PUV  #LE«  PUVMXi  GO  TO  10 

V0bU0027 

WFITf  (  6*90)  PUV 

VOwUOOSO 

GO  TO  18 

V0vU0033 

10 

veo  «  puv  ♦  9o«oo 

VOwUOOf* 

REO  ■  A.O  ♦  .04  «  TUFLW 

V0kU0039 

IF  (  RED  »LE*  REOMXI  GO  TQ  19 

V0wU0042 

WRITE  (6*91)  REO 

V0wUC045 

GO  TO  18 

VOkUOOAf 

19 

REL  ■  9.0  ♦  .10  •  TUFLW 

VOwUOOSl 

IF  I  REL  «LE.  RELMX)  GO  TO  20 

V0wU0094 

WRITE  (6*921  REL 

VOwUOOS? 

10 

IERFLG  •  1 

VOkUOOBO 

TOTV  >  0.0 

VOi.U0063 

RETURN 

V0*.U0066 

20 

ACTL  ■  2Z1  *  R  4  13.0 

V0kU0069 

IF  (  IFPFL6  .NE.  01  GO  TO  21 

V0wU0072 

VEQN  ■  2.C  *  VEO 

V0-U0079 

REON  «  2.0  •  RED 

V0bU0078 

RELN  -  2.0  •  REL 

VOwUOOBl 

ACTLN  ■  2.0  *  ACTL 

VOwUOOBA 

21 

I F < 9.6  -  R 1  22*23*24 

22 

SGN  •  -1. 

GO  TO  29 

23 

SGN  *  0. 

GO  TO  29 

24 

SGN  «  1. 

29 

TOTV  ■  VEO/VEON  ♦  RED/REDN  ♦  REL /RELN  ♦  ACTL/ACTLN  ♦ 

*123«(1.-SGN» 

IF  (IFPFLG  .EO.  0)  TNOMV  •  TOTV 

VOuUOOS* 

RETURN 

VOvUOOSO 

90 

FORMAT  UH  «  10X*  6HPUV  >  F8.2*  12H  EXCEEDS  MAX  ) 

VOuUOOfl 

51 

FORMAT  (1H  *  10X*  6HRED  -  F8.2*  12H  EXCEEDS  MAX  > 

V0bU0096 

92 

FORMAT  (1H  »  10X*  6HREL  ■  F8.2.  12H  EXCEEDS  MAX  ) 

END 

VOwU0099 
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-IT  FOR  RELIAB 

subroutine  RELIAB 

COMMON  /DATA/  DUMU9S.  RMAX 
COMMON  /FIG/  IFPFLG.  f ERFLG 
COMMON  / I NTM /  DUMl(4»t  0.  PRES 
COMMON  /OUR/  DUM2<3>.  TOTR 
COMMON  /NOMI/  DUM3 { 3 ) .TNOMR 
IF  ( IFPFLG  • NE i  0)  GO  TO  20 
ACCR  «  12»3 
RER  *  16.0 
TUR  «  44.0 
VAPR  «  16.0 

PTOTR  «  ACCR  +  RER  *  TUR  ♦  VAPR 
20  PUR  *  463.0  ♦  0.0111  *  PRES  *  0 

TOTR  «  PTOTR  ♦  PUR 
IF  (  IFPFLG  .FQ.  OS  TNOMR  >  TOTR 
IF  (  TOTR  .IE.  RMAX)  GO  TO  30 
IFRFLG  «  1 
WRITE  (6.40) 

30  RETURN 
40  FORMAT 
ENO 


REkI 0000 
REkt  0003 
REkt  0006 
REu 1 0009 
REbt00l2 
RE*,!  0015 
R£kI0018 
REwIOOjl 
REkI 0024 
REwI0027 
REkI 0030 
RFwI0033 
REwI0036 
REwI0039 
REk 10042 
REk!0043 
REktOCAS 
REkt  0051 
REkI0054 
REwI 0057 
REuI 0060 


C 1H  .10X.23HR ELI  ABILITY  EXCEEDS  MAX  S 


XQT  MAIN 


5000. 0 

2000.0 

7.5 

250.0 

“0.5 

1.0 

0.75 

1.0 

0.75 

0.75 

0.625 

0.75 

0*625 

99999.999 

0.0 

90999.999 

0.0 

99999.999 

99999.999 

99999.999 

0.0 

40.0 

0.0 

0.0 

0.0 

185.0 

40.0 

185.0 

40.0 

5.0 

4C.0 

0.0 

48000.0 

99999.999 

•30.0 

15.0 

25.0 

25.0 

dataoooi 

OA  TA0002 
DATAOOOI 
DATA0004 
DATAPC05 
0ATA0006 
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APPENDIX  IV 

NO.  2  COMPUTER  PROGRAM 


1.  RANDOM  AND  ADAPTIVE  RANDOM  SEARCH 

Before  proceeding  to  the  details  of  the  computer  mechanization,  a  brief 
mathematical  description  of  the  searching  techniques  is  given.  More  details 
ere  presented  in  Reference  8. 

2.  RANDOM  OPTIMIZATION  METHOD 

The  random  optimization  method  will  be  described,  using  tl  e  notation  in  the 
statement  of  the  problem  in  Section  III,  If  u(l)  denotes  the  initial  parameter 
set  of  the  system,  the  performance  index  corresponding  to  this  state  will  be 
denoted  by  P,  I.  (1).  Let  £  be  an  n-dimensional  normal  random  vector  with 
zero  mean  value  and  unit  correlation  matrix.  Compute  the  P.  L  at  the  ran¬ 
dom  parameter  vector  u(l)  +  £(1)  and  call  the  corresponding  random  step  a 
success  if 


P.I.  [u(i)  +£(l)]  >  P.L(i)  +«. 

for  an  *_  chosen  a  priori.  Otherwise,  the  random  step  is  called  a  failure. 
Depending  on  whether  the  first  random  step  is  a  success  or  failure,  the 
parameter  vector  u(2)  is  defined  as  follows: 

u(2)  =  u(l)  +£(1)  fora  success 

u(2)  s  u(l)  for  a  failure 


Next,  the  performance  criterion  is  computed  at  the  random  parameter 
vector  u(2)  +  £(2),  and  the  procedure  is  repeated.  The  recursive  relation¬ 
ship  atthe  kfh  random  step  is  seen  to  be:  (1)  compute  P.  L  [u(k)  +  £(k)], 

(2)  is  P.  L  [u(k)  +  4(h)]  >  P.  L  [u(k)]  +  »_,  (3)  if  step  2  is  a  success,  define 

u(k  +  1)  =  ufr)  +  £(k),  if  not  u(£  +  1)  =  u(k).  By  continuing  this  process,  a 
random  sequence  of  states  u(k)  is  obtained.  The  convergence  of  this 
sequence  to  the  optimal  state  is  guaranteed. 

3.  ADAPTIVE  RANDOM  OPTIMIZATION  METHOD 

The  adaptive  random  optimization  method  is  essentially  the  same  as  the 
method  just  described.  To  speed  the  convergence  to  the  optimum  system, 
however,  the  random  vector  is  chosen  in  a  slightly  different  manner.  The 
mean  of  the  random  vector  and  its  variance  will  be  varied  according  tc  past 

a  MM  /«  a  /%  f  sVa  amS!  n  •  f  I  m«i  Am  a  ••  S  tm  iJafi  •*  •  *  •  m  <4  i*>  A  m  a  ■  a 

V  A  *VHVV  VS  V**  V  vj#«at*4SSiiS«SV(4  V  f  v  •  *  (  WVSSMV  w  *  »•»  |SS  W  • 

£{k)  o  d(k)  +  T(k)i(k) 

where  d(k)  is  the  mean  value  of  t(k),  and  T(k)  is  a  transformation  matrix 
operatfng  on  the  random  vector  J(k). 
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Once  the  random  process  is  defined,  it  is  used  in  the  same  manner  as 
the  random  vector  was  used  in  the  random  search  method.  All  that  remains 
is  todescribe  calculation  of  the  random  process.  If  the  first  step  is  to  be  taken 
in  a  truly  random  manner,  the  initial  value  of  the  mean  value  vector  should 
be  zero.  If  the  initial  step  was  a  succees,  then  the  mean  value  of  the  second 
step  is  chosen  to  bias  the  choice  of  the  random  step  in  the  same  direction. 

That  is, 


d(Z)  =  Cjid). 


If  the  initial  step  was  a  failure,  the  mean  value  of  the  second  step  remains 
zero  or  is  chosen  so  as  to  ensure  motion  in  the  opposite  direction.  The 
recursion  relation  for  the  k^1  random  step  is  defined  to  be 


d(k  +  1)  =  cQd(k)  +  c,£(k) 

where  cn  and  c,  satisfy  the  following  conditions 

v  a 


CSc0<<,  c,!0,  c0  +<:,>! 


when 

PI(k)  >PI(k  -  1)  +  f_ 

0  <  CqC  1.  Cj  <  0,  /cQ  +Cj/  <  1 

when 


PI(k)  <  PI(k  -  1)  +  f. 


Updating  the  mean  value  in  the  manner  described  above  increases  the  prob¬ 
ability  of  a  successful  step,  because  it  coincides  vnth  the  direction  of  expected 
increase  of  the  performance  index.  It  also  lengthens  the  random  step  and 
speeds  up  the  process  of  finding  the  optimal  system.  Douglas  has  experienced 
success  in  applying  this  technique  and  has  developed  an  algorithm  for  updating 
variance  of  the  random  process  similar  to  the  technique  described  above  for 
updating  the  mean  value.  The  results  of  this  work  are  discussed  in 
References  )  and  Z. 

4.  DESCRIPTION  OF  THE  ADAPTIVE  RANDOM  SEARCH  ROUTINE 


The  cumyutci  prGgiiiT*  used  essentially  the  same  format  that  was  used  in  the 
fixed-grid  program.  The  subroutines  used  to  evaluate  the  cost  functions 
were  identical.  The  main  program  was  changed  to  incorporate  the  random 
search  technique.  In  addition,  it  was  found  to  be  simpler  to  incorporate  the 
design  of  the  nominal  system  into  the  main  program. 


The  flow  diagram  chown  in  Figure  13  shows  the  logic  of  the  routine.  The 
listing  that  follows  (No.  2  Program)  is  only  for  the  main  program  and  shows 
how  the  mathematical  description  of  the  adaptive  random  search  routine  was 
incorporated  into  the  program.  Random  numbers  were  selected  for  all  of  the 
independent  variables  (pressure,  area,  moment  arm,  and  four  tube  sices). 
Given  this  selection,  this  set  had  to  be  tested  to  ensure  that  it  was  compatible 
with  the  design  requirements,  such  as,  static  hinge  moment,  flow  at  -40*F, 
and  dynamic  hinge  moment.  If  all  of  the  tests  were  passed,  the  parameter  set 
was  used  to  calculate  the  P,  I.  ;  if  not,  a  new  eat  of  parameters  were  selected 
and  the  process  repeated.  The  old  set  was  classed  as  unsuccessful  for  the 
purpose  of  updating  the  adaptive  feature  of  the  routine. 

Unsuccessful  trials  were  counted  at  three  places  in  the  program:  (1)  if  the 
hinge-moment  test  (ARTEST)  failed;  (2)  if  the  -40*  F  flow  requirement  'ailed; 
(3)  if  the  P.  I.  failed  to  exceed  the  stored  P.  I.  by  *.  If  the  count  in  ar.)  of 
these  threa  places  exceeded  constants  chosen  a  priori,  the  run  was  halted. 
Tests  1  and  2  were  incorporated  to  halt  the  program  If  it  reached  a  situation 
where  it  could  not  pass  any  of  the  tests.  After  the  initial  program  was  oper¬ 
ating  properly,  this  situation  was  never  encountered. 
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LISTING  FOR  NO.  2  COMPUTER  PROGRAM 


OPTIMAL  0ES1GN  OP  AN  ACTUATING  SVSTEM  KAIN0000 

DIMENSION  FJSU|,RW<7l,0L00m,0LDXm,0(7l 

COMMON  /OATA/  AI381  MAINS006 

COMMON  /NOMI/  TNOMW  »TNOMO»  TNOMV, TNQNR  *TNQMOV  • TNOMEN  MAIN9009 

COMMON  /OUP/  TQTW » T OTC » TOT V, TU TR *  TOT DVP, TOTE NV  MAIN00I2 

COMMON  /FIG /  IFPFLG, IERFLG 

COMMON  /INTrt/  DSl  t  OS2  >OAL  1  •  DAL  2.  Q,  PRFS,  ARE A, R , ACTOP , TUFLW, ANORM  MAIN001S 
EQUIVALENCE  ( At  2  I «PM) , ( A( 41 , P8» , I A(  5 1  ,RN» , 

1(  At  15J,A6*,(A(  l 7), AMI. (A( 231  , ABI ,  I  A(  26 1 ,TULAL 1) , ( A< 27 1 , TULAL2I , 

2t  A(  28} » TUI  5 1 1 ,( AI29 1 ,  TULS2 I , ( A ( 30) «  DEL  MX  I  ,(A(3ll,HM> 


WRITE  (6,553)  MMN0023 
READ  (5,5001  (All), l-l, 38*  MAIN0024 
WRITE  16,550)  (All), 1-1,381  MAIN0027 
WRITE  (6,5551  MAIN0034 
IFPFLG  -  0  MAIN0042 
IERFLG  •  0  MAIN0045 


NARK  •  15 
PIi)i  •  0.0 
NO  -  0 
NUN  -  0 
KOUNT  •  0 
PRES  •  A(  II 
R  -  A  (  3  I 
DSl-A(8l 
0S2-AI 121 
OAL l»A( 6) 
0AL2-AI10) 


AREA  • 

MM/I  PRES*  A  1 

NOMS0021 

ANORM 

•  AREA 

NOMS0121 

228  - 

PRES  -  3000. 

IFI228 

.LT.O.OI  228  -0.0 

221  - 

.001705*1  TULS1*(0SI>*«(- i .931*11. *.0002*2281**2  ♦ 

TULAL1*( 

10ALI ‘**(-6.9311 

222  • 

.001 705*<  TULS2*( DS2 !*•( -4.931*1 l 0002*2281 **2  ♦ 

TULAL2* I 

10AL2I**«-4.931 1 

225  - 

•004533*0ELMX 

NOMSOOll 

Q  *  225*AREA«R 

N0MS00V3 

223  • 

14.4*  Ql ** 1 . 75 

NOMS0045 

TUOP 

•  221*18.8*01**1.75  *  222*223 

NOMS0048 

TU0P2 

-  221*  223  ♦  222  *(2.2*QI**I.75 

NOMS0051 

ACTOP 

•  PRES  -  .25*1  PRES  -  TUOP  1  -  IU0P2 

N0MS0054 

224  • 

I  .75*PRES  -  ACTOP 1 

IF  1 224. LE. 100.01  GO  TO  30 

NQMS0Q6/ 

PRES  - 

PRES  *  224/  2.0 

NOMS0063 

GO  TO 

2 

NOMS0066 

00  16 

J-1,7 

SUM  12  -  0.0 

00  15  (-1,12 

CALL  RANOMl  I  MARK, X) 

RANO  •  X 

SUM  12  •  SUM12  ♦  RAND 
!«,  CONTINUE- 

RN( J I  •  SUM 12  >6.0 

0(J1  •  CO*  WOO(J)  ♦  Cl  *  OLDX(J) 

OLDOU)  «  OUI 

RN(  J I  -  OUI  *  RN(J) 

QLDXU)  -  RNU! 

16  CONTINUE 

PRES  -  PRES  »  RN( 1 1  *  PM 
R  •  R  *  RN12I  *  RM 
AREA  •  AREA  «-  RN(3)«  AM 


ARTEST»PRES  *Ar<EA*R 
IFIARTEST.LT.HM)  CO  TO  220 
FART  •  0.0 
OSl-ACGi 

IFIRNCAl.LT. PARTI  DSI-AC9} 

DS2-AC 121 

IFCRNC9I.lt. PARTI  0S2>AC13I 
0AL1-AC6) 

IFIRNC6I.Lr. PART)  OAU-ACT) 

0AL2-ACI0I 

CFCRNC7I.lt. PART)  0AL2-AIU) 

23  ZZ1  -  .3366*1  ?ULSI*(0S1 )**( -3 .981 *( l .*.0002*2281**2  *  TULAL1*C0AL 
1  i I •• C  —  3 .9S i I 

112  -  . 3366*C  TULS2*C0S21**C -3.981*1 1.* .0002* ZZ 81**2  *  TULAL2*C0AL 
12I**C-3.98»1 

ZZ9  »  .Q0A9J3*0£LMX  NQMSOOll 

Q  •  ZZ9*AREA*R  NOMSOOA3 

29  7U0PA  «  Q*  f  2.2*111  ♦  ZZ2I  NAIN01S6 

(FIPRES.LT. TUOPAI  CO  TO  222 

ZZI  -  .00I709*(  TULS1*C  0$11**( -A .931*11. *.0003*Z28)**2  *  TULAL1*( 

10ALI l**(~A.931 I 

112  •  .001709*1  TULS2*(DS2)**C-A.93)*( l.*.0003*ZZ8l**2  *  TULAL2X 

!0 AL  2  )**(*> A.  93 1 1 
ZZ3  -  (A. A*  Q) **l. 79 

TUOP  -  Z 2 1 *( 8 .8*0 1 **1.79  ♦  ZZ2*2Z3 
TU0P2  ■  ZZI*  ZZ3  ♦  112  *C2.2*0l**i.75 
ACTOP  •  PRES  -  .25*C  PRES  -  TUOP  I  -  TUDP2 
CFC  C .7S*HM-1000.1-ACT0P*AREA*A1  30.30,396 

30  WRITE  C6.99AI  PRES, AREA.  Q,  TU0P.TU0P2, ACTOP, R 
WRITE  16,960)  COCCI,  I*>1,  Tl ,  CRN!  Jl,  J«1  *  71 
NO  •  0 
NUN  «  0 
CALL  WEIGHT 

CFC  lERFLG  .EQ.  01  GO  TO  39 
32  PIOX  «  0.0 

lERFLG  -  0 

IF  (IFPFLG  .EQ.  0)  GO  TO  210 
GO  TO  396 
39  CALL  COST 

IF  C  lERFLG  .NE.  0)  GO  TO  32 
CALL  OVNPER 

IF  I  lERFLG  .NE.  01  GO  TO  32 
CALL  VOLUME 

IF  C  lERFLG  .NE.  0)  GO  TO  32 
CALL  RELIAB 

IF  (  lERFLG  .NE.  0  I  GO  TO  32 
IFPFLG  »  IFPFLG  ♦  l 


110  PIOX  -  0.0  MAIN0193 

F JSC  1 1  «  1.0  -  TQTW/C 2»0*TNONH)  NA'NOISA 

F JSC 21  ■  1.0  -  TOTC/C  2.0*TN0WC I  MINI  ISA 

FJSC31  •  l.  -  TOTOYP 
FJS(**lm  l.O— T0TV/(2->0*TN0NVI 
FJSC9)  •  1.0  -  TOTR/  (2.0  *  T NO MR  I 

00  112  II  •  1,9  NAINA ISA 

IF  (CFJSC III.LT.  O.Ol.OR.CFJSCIII.GT.l.OII  GO  TO  119  KAIN7194 

*  *2  CONTINUE  NAINA ISA 

00  113  II  •  1,9  MAIN919A 

113  PIOX  -  PIOX  ♦  AC  11*31 1 *F JSC  III  NAIN0199 

POT*  PIOI  *  0.02 
IFCPIOX.GF.POTI  GO  TO  UA 
CO  >0.6 


NOHSOOAS 

NQNS00A8 

N0NS0CS1 

NONSOOSA 


NA3N0138 

NAIN01A1 

NAIN01AA 

NAIN01A? 

NAIN01A8 

NA I NO 190 
NAIN0192 
NAIN1192 
NAIN2152 


111 


Cl  *  -0.7 

KOUNT  »  KOUNT  *  1 

PRES  -  PRES  -  RN(l)  «  PN 

R  -  ?  -  RN( 21  *  RN 

AREA  «  AREA  -  RN(3)  *  AN 

I F( KOUNT. EC.  500)  GO  TO  210 

GO  TO  115 

114  PIOI  >  PIOX 
CO  -0.9 

C)  *0.2 
KOUNT  *  0 

115  11*1 

(FIOSl.EQ. 0.750)  11*2 
I  J*1 

1E10S2.EQ. 0.6251  IJ-2 
IK«1 

IFIDALl.EQ. 0.750)  IK*2 
IL-1 

(El DAL2.EQ.O.tj59  IL*2 

WRITE  16,552)  U,U,!I,IK,  PIOX  ,  TOTW,  TOTC,  TOTOYP,  TOTV,  TOTR 
IF  I IFPEIG.NE. 1)  GO  TO  4 
PIDI-0.0 
00  165  K  •  1,7 
OLOOIK)  -  0.0 
OLOXIK )  -  0.0 
165  CONTINUE 

PRES  •  6917. 

AREA  *  3,19 
R  -  7.7 
GO  TO  23 

WRITE  16*551)  PIOI 
STOP 

220  WRITE  (6,556)  PRES, AREA, R 
NO  -NO*- 1 

WRITE  (6,560)  ( 0(  1 1 , 1*1, 7) , IRN( J ) , J-l , 7) 

IF( N0.NE.60l  GO  TO  356 

221  STOP 

2 22  (1-1 

IFI0S1.EQ. 0.750)  1 1 -2 
I  J-l 

IF( 0S2. EC. 0.625)  IJ-2 
IK-1 

IFlDALl.EQ. 0.750)  IK-2 
IL-l 

IF(0AL2.EQ. 0.625)  IL*2 

WRITE  (6,5571  IJ,1L,II,IK,  PREf.,AREA,R 

NUN  -NUN  *  1 

IFINUN.NE.40)  GO  TO  956 

GO  TO  221 

356  PRES  *  PRES  -  RNI1I  *  PN 
R  »  R  -  RN(2)  •  RN 
AREA  -  AREA  -  RN(3)  *  AN 
CO  -0.6 
Cl  •  -0.7 
GO  TO  4 

«'*')  FORWAT  I  7F10.3  ) 

J  FORMAT! 1H0.47X, 10H1NPUT  DATA////13H  P  NON  •  F10.3,«X  PN  • 
IF  10. 9, 10H  R  NON  -  F10.3.8H  P8  *  F10.3,3H  RN  • 

2F1G.1//13H  0AL11  •  F10.3,11H  QAL12  •  F10.3,10H  0S11  • 

3F  10*3, 10H  0S12  -  F1Q.3,UH  0AL21  •  F10.3,11H  0AL22  *  F10. 

413H  0S21  •  F10.3.UM  OS22  •  F10.3.10H  W  MAX  •  FlO.?  , 


NAIN01S>> 


NAIN0195 


NAIN0199 
NAIN0204 
NAIN0207 
NAIN022G 
3//HAIN02 13 
NAIN0216 
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58H  AB  -  F10.3.UH  C  MAX  »  FI0<,3,JH  AN  •  *10.3,  //13H  P'JV  NAXMAIN0219 
6  <  F10.3, IIH  REO  MAX  -  ?\0.3,lOH  «  MAX  •  FI0.3.10H  R  MIN  »  NA 1*0222 

7F10.3*11H  OYP  MAX  »  FIO.3,  8H  AB  -  FIO. 3,  //13H  ENV  MAX  •  MAIN0225 

8F10.3.11H  ENV  MIN  •  F10.3,  OH  TULAL1  «  F10.3*10H  TULAL2  *  FIO.3,  MAIN0228 

9UH  T'JtSl  »  F 10* 3#  1 IH  TUIS2  >  F10.?//13H  DELTA  MAX  •  F|0.3t  NAIN0231 

lllrt  HM  »  FIO.3, IOH  WTFk  -  FI0.3, 10H  C0STF2  -  FI0.3, IIH  0YPFMAIN0232 

23  -  FIO.  3//13H  VOLF  A  -  .10. 3, IIH  RELf-5  «  F10.3,10H  ENVF6  NAIN1232 

3>  F 10.3, 1 IH  REL  MAX  •  FIO. 3)  MAIM2232 

551  FORMAT! 1HC, ATX, 28H  THE  STORED  VALUE  OF  P.l.  •  FIO. 31 

552  FORMAT  (114  ,  AH  CASE  AI  1,6X,  7HP.I.  *  F7.2,3X,5HVY  »  F8. 1 ,3X,  TMCJST 
I  «  F7.0,  3X,6H0YP  •  F6.2 , 3X, 6HV0L  -  F6.2,3X,6K*EL  •  FT. 1 i 

553  FORMAT  (1HI,  AOX,  2AH0PT INAL  ACTUATING  SYSTEM  I  MA1N02A1 

55A  FORMATUHO,  8H  PRES  »  FIO.  3,  8H  AREA  «  FlvJ.3,  5H  Q  «  FIO.3,  NOKS0075 

19H  TUOP  ■  F 10. 3, 1  OH  TU0P2  >  F10.3.9H  AC  TOP  •  F10.3,5H  R  FIO.3! 

555  FORMAT! 1HI,3AX,25HCASE  »  0S2,0AL2,DSI ,DAL1  I  NAINI2A3 

556  FORMAT  I1H0,8H  PRES  »  FIO. 3,  8H  AREA  -  FIO.3,  5H  R  *  FIO. 31 

557  FORMAT  (1H0,6H  CASE  AI1,6X,8H  PRES  -  FIO.3,  8H  AREA  -  FIO.3, 

15H  R  »  FIG. 31 

560  FORMAT  (IH  ,3H7US  7F8.2.5X.AH7RNS  7F8.2I 

ENO  MAIN02AA 


133 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  RAD 

(Imaurlty  eiftaaitieotion  o  t  titia.  body  ol  abotroct  our*  fndotdng  annotation  mm  at  bo  anttrod  w*on  N  ooaroli  rtjMif  io  cl aaaitiod) 


1 .  OM'  KAYIN  O  ACTIV|TY  (Corporota  author)  Sc.  IIRORT  SECURITY  C  LASSiriCA YION 

Douglas  Aircraft  Company,  Inc.  Unclassified 

Missile  and  Space  Systems  Division  it  anoua - 

Santa  Monica,  California  N/A 


1.  RERORT  TITLE 

Optimal  Actuation  Research  and  Study 


4.  DCSCniPTIVE  NOTES  (Typo  of  report  and  inctuaira  data*) 

Final  Report  -  1  March  1966-1  February  1967 


S-  i^ltTNONf ;Lm(  ncmc.  tint  nmmo.  initial) 


Price,  Harold  A. 
Giesekine,  Darrell  L. 


«.  REPORT  DATE 

.tuly  1967 


it  CONTRACT  OR  SNANT  NO 

AF33(615)-3514 


7*.  total  no.  of  pami  7*.  no.  or  rin 

|  140  B 


»•.  oRiaiNATOR’t  report  NuuaeRfSj 


8225 

c  Task  No. 
822510 


DAC  -  60576 


*k  JThim^poat  NOTJ;  (Anythin 

AiFDL-TR-67-46 


i  An  mmy  a*  miaiN 


10  A  VAIL  AEILIVY/LIMITATION  NOTICE* 

This  document  is  subject  to  special  export  controls  and  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with  prior  approval 
of  the  Air  Force  Flight  Dvnamics  Laboratory,  Air  Force  Systems  Command, 

1 iVl ^Vl4i&iMrY  N&li*  f°rce'  It  IPONIORINO  MILITARY  ACTIVITY 

.  Air  Force  Flight  Dynamics  Laboratory 

'  Wright-Patterson  AFB,  Ohio  45433 


1».  ABSTRACT 

The  servo  actuating  and  corresponding  energy  distribution  subsystem  com¬ 
prise  approximately  80%  of  an  aircraft's  flight  control  system  weight.  Conse¬ 
quently,  whenever  subsystem  tradeoff  studies  are  conducted  it  is  desirable  to 
select  the  optimum  design  with  respect  to  weight  and  other  similar  parameters. 
This  study  investigates  and  develops  such  an  optimal  design  process.  A  sample 
problem  was  selected  and  an  optimal  technique  formulated  and  applied  to  the  i 
problem.  The  sample  problem  was  a  fixed -configuration  hydraulic  actuation  and  I 
power  system.  The  study  objectives  were  to  optimize  weight,  dollar  cost,  size, 
dynamic  performance,  and  reliability  as  a  function  of  the  system's  independent 
design  parameters.  The  parameters  included  pressure,  actuator  area, 
actuator  torque  arm,  and  plumbing  tube  sizes.  Parameter  optimization  was 
accomplished  by  fixed  grid  and  random  searching  techniques.  Within  the  frame¬ 
work  of  parameter  optimization,  a  design  philosophy  was  formulated  which 
allowed  dissimilar  terms  (e.  g.  ,  weight  in  pounds  and  dollar  cost  in  dollars)  to 
be  combined  io  form  a  xoxai  performance  criterion  for  the  system.  When  the 
optimization  technique  was  applied  to  the  sample  problem,  the  performance 
criterion  showed  little  variation  as  a  function  of  the  parameters  being  optimized. 
All  of  the  cost  functions  had  large  nominal  values  and  only  slight  variations 
about  this  nominal.  To  realize  the  potential  of  the  design  technique  developed 
in  this  study,  different  design  concepts  and  possible  variations  of  each  one  should 
be  considered,  to  reach  a  more  meaningful  optimum  design. 


DD  MX.  1473 


Security  CUsslflcstioa 


Unclassified 


Security  Classification 


KEY  *0*0* 


1.  Actuation  System 

2.  Flight  Control  Actuation  System 

3.  Optimal  Design 

4.  Optimal  Actuation  Systems 


LINK  A 


LINK  B 


LINK  C 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Em«r  the  nema  and  aMrtii 
of  tha  contractor,  aubcontractor,  grantaa,  Department  of  Da- 
fanaa  activity  or  other  organization  (corporal*  author)  issuing 
the  report. 

2a.  REPORT  SECUHTY  CLASSIFICATION:  Enter  the  over¬ 
all  aacurity  claarification  of  the  report.  Indicate  whether 
"Reatricted  Data"  ia  included.  Marking  la  to  be  in  accord¬ 
ance  with  appropriate  aacurity  regulation*. 

2b  GROUP:  Automatic  downgrading  ia  opacified  in  LoD  DP 
rectiva  5200. 10  and  Armed  Forcaa  Industrial  Manual.  Enter 
the  group  number.  Aleo.  when  applicable,  ahow  that  optional 
markings  have  been  ueed  for  Group  3  and  Group  4  a*  author¬ 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  let'em.  TIUea  in  all  caaea  ahottld  ba  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  ahow  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES  If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  summery,  annual,  or  final. 

Give,  the  incluaiva  dates  when  a  specific  reporting  period  la 
covered. 

5.  AUTHORISE  Enter  the  name(s)  of  author!  a)  as  a  ho  am  on 
or  in  the  report.  Enter  last  name,  flrat  name,  middle  initial. 

If  military,  ahow  rank  and  branch  of  service.  The  name  of 
the  principel  author  ia  an  abaoluta  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  data  of  the  report  at  day, 
month,  year,  or  month,  ysar.  If  more  than  ona  data  appears 
on  the  report,  use  deta  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures.  La.,  enter  the 
number  of  pages  containing  information. 

7u.  NUMBER  OF  REFERENCES  Enter  the  total  number  of 
references  cited  in  the  report. 

8*  CONTRACT  OR  GRANT  NUMBER:  It  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

bb.  Sc,  &  8 d.  PROJECT  NUMBER  Entar  the  appropriate 
military  department  identification,  such  aa  project  number, 
subproject  number,  system  numbers,  task  number,  ate. 

9a.  ORIGINATOR'S  REPORT  NUMBER'S):  Eater  the  offi¬ 
cial  report  number  by  sehlch  the  document  will  be  Identified 
and  controlled  by  tha  originating  activity.  This  number  must 
ba  unique  to  this  report. 

96.  OTHER  REPORT  NUMBEJKS):  If  tha  report  has  been 
assigned  any  other  report  numbers  farther  by  th a  originator 
or  by  the  rponrjt),  also  enter  this  number! a). 

10.  AV AIL ABn .IT Y /LIMITATION  NOTICES  Enter  any  llm- 
nations  on  lurther  dissemination  of  the  report,  othar  than  those 


imposed  by  aacurity  classification,  using  standard  statements 
such  at: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC." 

(2)  "Foreign  announcement  end  dissemination  of  this 
report  by  DDC  ia  not  authorized  " 

(3)  ”U.  S.  Government  agencies  may  obtain  coplea  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
user*  shall  requaat  through 


(4)  "U.  S  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC.  Other  qualified  users 
shall  request  through 


(5)  “All  distribution  of  tnis  report  la  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 


If  the  report  has  been  furnished  to  the  Office  of  Technical 
Service*,  Department  of  Commerce,  for  aala  to  tha  public,  indi¬ 
cate  this  fact  and  entar  the  price,  If  known> 

11.  SUPPLEMENTARY  NOTES:  Gee  for  additional  explana¬ 
tory  notaa. 

11  SPONSu  inG  MILITARY  ACTIVITY:  Enter  tbs  name  of 
tha  departmental  project  office  or  laboratory  sponsoring  (pay¬ 
ing  lor)  the  research  and  development.  Include  address. 

13.  ABSTRACT:  Entar  an  abstract  giving  a  brief  and  factual 
summary  of  tha  document  Indicative  of  tha  report  even  though 
it  may  also  cprear  elsewhere  in  tha  body  of  tbs  technical  re¬ 
port.  If  additional  apaca  la  required,  a  continuation  sheet  shall 
ba  attached. 

It  is  highly  desirable  that  tha  abstract  of  classified  report* 
be  unclassified.  Each  paragraph  of  tha  abstract  shall  and  with 
an  Indication  of  tha  military  security  classification  of  the  in¬ 
formation  in  tho  paragraph,  represented  aa  ( Tt).  (t).  (C).  or  (V). 
Thin  la  no  limitation  on  tho  tmirHi  of  the  ahetreet.  Kow. 

avar,  tha  suggested  length  is  from  1$5  to  225  won]*. 

14.  KEY  VCRDS:  Kay  words  are  technically  meaningful  tanas 
or  short  phrasas  that  characterize  a  report  and  may  be  uied  ao 
lcdai  entries  for  cataloging  tha  report  Key  words  must  b* 

••lac tad  so  that  no  aacurity  classification  ia  required.  Identi¬ 
fiers,  such  as  equipment  mottal  designation,  trad*  name,  military 
project  coda  nama,  geographic  location,  may  be  used  as  key 
words  but  edil  be  followed  by  an  indication  of  technical  con- 
tait.  Tha  assignment  of  links,  rules,  and  weights  1*  optional. 


Unclassified _ 

Security  Classification 


